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Chapter  1 

DISCUSSION  OF  SIGNIFICANCE  AND  SCOPE 


OF  THE  UNDERTAKEN  STUDY 


1 .1  Introduction 

Research  and  development  in  advanced  composites  technology  and 
design  procedures  has  advanced  to  a  point  where  these  materials  can 
now  meet  the  challenge  of  highly  demanding  aerospace,  automobile  and 
other  structural  applications.  Performance  payoffs  in  the  form  of 
weight  savings,  component  integration  and  a  five-fold  maintenance  re¬ 
duction  as  compared  to  metal  design  have  been  demonstrated  (Forsch 
[1])  and  contribute  to  the  overall  picture  of  an  attractive  marketable 
product. 

However,  resin  matrix  composites  have  moduli  and  strength 
properties  which  are  viscoelastic  or  time  dependent.  As  a  result  there 
is  a  need  for  methods  by  which  short  term  test  results  can  be  used  to 
predict  long  term  results. 

The  current  effort  is  part  of  a  continuing  cooperative  NASA-Ames  - 
VPI&SU  research  program  directed  towards  the  development  of  an 
accelerated  strength  and  stiffness  characterization  procedure  which 
generates  the  design  data  and  the  confidence  levels,  needed  to  insure 
structural  integrity  for  any  projected  lifetime.  This  process  is 
greatly  complicated  by  a  number  of  environment  exposure  variables 
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such  as  temperature,  moisture*,  UV  and  space  radiation,  high  vacuum 
and  many  others.  An  example  of  a  typical  thermal  profile  for  an  air 
carrier  given  by  Ripley  [2]  is  given  in  Fig.  1.1.  Any  attempt  to 
handle  such  variables  by  overdesign  will  penalize  the  economical 
features  mentioned  above  as  a  knockdown  factor  will  seriously  limit 
the  full  weight  reduction  potential  of  advanced  composites. 

1 .2  Previous  Efforts 

An  accelerated  characterization,  which  is  schematically  given  in 
Fig.  1.2,  has  been  developed  and  documented  by  Brinson  et  al .  [4-13], 

A  fair  amount  of  experimental  support  has  been  accumulated,  based  on 
elements  of  linear  viscoelasticity  by  Yeow  [5].  Nonlinear  visco¬ 
elasticity  came  into  this  picture  when  Griffith  [7]  potentially  identi¬ 
fied  this  behavior  as  the  reason  for  a  serious  underprediction  of 
creep  rupture  data  for  [90/±60/90]2s  laminates.  The  prediction  was 
based  on  a  linear  incremental  lamination  theory  developed  by  Yeow  [6], 
As  a  result,  a  study  of  nonlinear  viscoelastic  characterization  pro¬ 
cedures  was  performed  by  Griffith  [7]  which  led  towards  a  valid 
graphical  Time-Temperature-Stress  Superposition  (TTSSP)  principle  for 
the  generation  of  modulus  master  curves. 

Dillard  [9]  developed  a  non-linear  viscoelastic  incremental 
lamination  theory  and  found  reasonable  correlation  between  his  predic¬ 
tions  and  experiments.  Because  of  the  limited  experimental  data 


*It  has  been  shown  that  moisture  pickup  is  generally  reversible. 
When  combined  with  thermal  spikes,  however,  this  phenomenon  can  create 
irreversible  damage  as  shown  by  Adamson  [3], 


Fig.  1.2  Accelerated  Characterization  Method  for 
Laminated  Composite  Materials. 
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available  from  the  work  of  Griffith,  Dillard  used  a  semi empirical  non¬ 
linear  viscoelastic  model  due  to  Findley  to  represent  lamina 
compliances. 

1 .3  Objectives  and  Outline  of  the  Present  Study 

The  objective  of  the  current  study  was  to  cast  the  graphical 
TTSSP  method  into  an  analytical  form.  It  was  believed  that  the 
analytical  framework  developed  by  Schapery  in  a  series  of  publications 
could  be  adapted  for  this  purpose.  A  literature  review  on  nonlinear 
viscoelastic  constitutive  equations  is  given  in  Chapter  2,  while  the 
thermodynamic  basis  and  the  derivation  of  the  Schapery  model  are  the 
subjects  of  Chapter  3. 

An  effort  to  relate  viscoelastic  deformation  and  failure  on  an 
energy  basis  is  developed  in  Chapter  4.  This  allows  an  immediate  trans¬ 
fer  of  any  modulus  master  curve  into  a  strength  master  curve  as 
depi cted  in  Fig.  1.2. 

Chapter  5  concentrates  solely  on  934  neat  resin  properties  and 
on  the  discussion  of  obtained  experimental  results  thereto. 

Master  curves  for  transverse  and  shear  properties  of  T300/934 
as  obtained  through  the  Schapery  analysis  are  discussed  in  Chapter  6. 

Results  on  creep  rupture  for  [90]g  and  [60]qs  lamina,  based  on 
the  energy  criterion  (Chapter  4),  together  with  a  comparison  of  an 
incremental  lamination  theory  prediction  with  experimental  results  on 
[90/±45/90]2s  are  discussed  in  Chapter  7. 

A  general  discussion  and  conclusions  are  given  in  Chapter  8. 
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Chapter  2 

INELASTIC  TIME-DEPENDENT  MATERIAL  BEHAVIOR 

2 . 1  Introduction 

The  purpose  of  this  chapter  is  to  produce  a  brief  literature 
overview  on  research  in  nonlinear  viscoelasticity.  The  majority  of 
publications  explain  how  to  deal  with  magnitude  nonlinearities,  i.e. 
they  account  for  the  stress  deoendence  observed  on  creep  compliance 
versus  time  as  illustrated  in  Fig.  2.1. 

A  much  smaller  number  of  articles  also  account  for  i ntermode 
non! i neari ty  (nonlinear  interaction  between  different  components  of  the 
stress  tensor)  and  almost  invariably  use  a  single  mechanical  in¬ 
variant,  namely  the  octahedral  stress.  In  this  respect  it  should  be 
mentioned  that  data  generated  by  Cole  and  Pipes  [14]  shown  in  Fig.  2.2 
indicate  the  extent  of  stress  interaction  in  off-axis  unidirectional 
boron-epoxy  coupons.  These  experimental  results  show  that  inelastic 
behavior  is  not  a  monotonic  function  of  fiber  orientation.  Similar 
behavior  has  been  observed  for  graphi te-polyimide  off-axis  coupons  by 
Pindera  and  Herakovich  [15],  They  indicate  that  the  inclusion  of 
residual  thermal  stresses  and  axial  stress  interaction  at  the  micro¬ 
mechanics  level  are  necessary  in  order  to  predict  the  correct  trends 
in  the  shear  response  of  off-axis  specimens. 

An  almost  never  mentioned  type  of  nonlinearity  is  interaction 
non! i neari ty ,  i.e.  nonlinear  interaction  between  events  (i.e.  load 
application  or  release)  occurrina  at  different  times. 
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Shear  Stress  ( Ksi  )  21  Stress 


g.  2.1  Stress  *  Strain  Behaviour  of  Elastic  and 
Viscoelastic  Materials  after  two  Values 
of  Elapsed  Time  ,  t . 


Shear  Strain 

Fig.  2.2  Extent  of  Stress  interaction  in  Off-axis 

Unidirectional  Boron-epoxy  Coupons- Cole  and 
Pipes  [14] ,  1974. 


Viscoelasticity 


2.2.1  Multiple  Integral  Representations 

The  multiple  integral  representation  is  a  very  appealing 
theoretical  concept,  since  it  is  not  limited  to  a  particular  material 
or  class  of  materials.  The  Volterra-Frechet  expansion  for  a  one¬ 
dimensional  case  is  written  as: 

rt  ft  ft 

e  =  0-j  (t  -  x i  ) da ( x i  )  +  D2(t  -  x-|  ,t  -  T2)da(x-|  )da(T2) 

j  —CO  '  -CO*'  —  oo 

+  ...  (2.1) 

The  first  integral  of  Eq.  (2.1)  describes  the  linear  viscoelastic 
behavior,  defined  by  the  Boltzmann  theory  and  the  second  and  higher 
order  integrals  are  representations  of  both  magnitude  nonlinearity  and 
interaction  nonl inearity,  the  latter  implying  an  interaction  effect 
between,  e.g.,  the  stress  increments  at  times  and  t2> 

If  a  force  is  suddenly  applied  to  a  specimen,  producing  a  constant 
stress  a,  it  follows  from  2.1  that 

e  =  D^tja  +  D2(t,t)a2  +  D3(t,t,t)a3  +  ...  (2.2) 

The  functions  D^(t  ,t , . . .  ,t)  can  thus  be  determined  from  creep  tests. 

The  number  of  stress  levels  must  be  equal  to  the  number  of  kernels  J^. 

Gradowczyk  [16]  showed,  however,  that  the  system  of  linear  equa¬ 
tions  which  serves  to  determine  the  unknown  functions  is  ill- 
conditioned  and  the  procedure  of  finding  the  kernels  is  unstable. 
Being  unaware  of  Gradowczyk 's  work,  this  author  attempted  to  use  the 
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multiple  integral  representation  for  p.v.c.  and  had  to  give  up  the 
attempts  due  to  the  extreme  instability  of  the  results  obtained  [19]. 

The  extension  of  the  Volterra-Frechet  expansion  to  the  3- 
dimensional  case  has  been  given  by  Green,  Rivlin  and  co-workers  [17]. 
Even  if  we  limit  their  expression  to  triple  integrals  (small  non- 
linearities)  twelve  kernels  still  must  be  determined  by  experiment. 

An  alternative  of  the  Green-Rivlin  approach  has  been  formulated 
by  Pipkin  and  Ropers  [18]  and  is  often  referred  to  as  nonlinear  super¬ 
position  theory  (NLST)  which  can  be  written  as 


-  t2  ,a(x2)  )ct(t-j  )ct(  t2  )dx-j  dt2  +  ...  (2.3) 

It  is  clear  that  eq.  2.3  is  nonlinear  even  on  the  first  approximation. 
Thus,  the  additivity  of  incremental  stress  effects  in  the  Boltzmann 
superposition  sense  is  preserved. 

Kinder  and  Sternstein  [20]  showed  that  difference  in  strains  pre¬ 
dicted  by  the  two  theories  is  in  the  transient  portion  if  the  material 
exhibits  fading  memory.  Thus, in  a  single  creep  test, the  equation  2.3 
reduces  to 

6(t)  =  (  D1 (t  -  T,3(T))5(-)dT  (2.4) 

J  -00 

which  is  the  most  general  representation  of  nonlinear  creep  provided 
the  transient  response  from  the  step  loadings  has  decayed.  However, 
whether  multiple  loading  and/or  unloading  programs  require  additional 
integrals  to  describe  time  interactions  of  the  loading  history  remains 
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an  open  question. 


2.2.2  Single  Integral  Representations 

i)  Modifications  of  the  Boltzman  Theory  of  Linear  Heredity 

Two  schools  of  thought  have  been  identified  in  the  literature. 
The  first  one  states  that  it  is,  in  principle,  possible  to  map  a  set 
of  nonlinear  isochronous  stress-strai n  curves  into  a  linear  set,  pro¬ 
vided  we  non! inearize  the  stress  and  strain  measures  as 


The  second  school  allows  for  a  stress  or  pressure-deformabl e  time 
measure.  Indeed,  the  pressure  dependence  of  the  mechanical  response  has 
been  predicted  theoreti cal ly  on  the  basis  of  the  free  volume  approach 
by  Ferry  and  Stratton  [24]  who  derived  the  pressure  analog  of  the  well- 
known  WLF  equation  in  the  form: 


log  ap  = 


(8/2.303  f0)(p  -  p0) 

VKf  -  -  p0J 


(2.8) 
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where  ap  is  the  ratio  of  the  relaxation  times  at  pressure  p  to  the 
relaxation  times  at  pressure  pQ,  B  is  a  constant  assumed  to  be  of  the 
order  of  unity,  f  is  the  fractional  free  volume  at  the  reference 
pressure,  and  is  the  isothermal  compressibility  of  the  free  volume. 
Beuche  [25]  and  O'Reilly  [26]  proposed  the  equation 

log  ap  =  C(p  -  pQ)  (2.9) 

where  C  is  a  constant  and  other  variables  are  as  defined  above. 

Knauss  and  Emri  [27]  postulate  a  linear  dependence  of  fractional 
free  volume,  on  the  instantaneous  values  of  temperature,  T;  water 
concentration,  C;  and  mechanical  dilatation  a,  i.e. 

f  =  f  +  aAT  +  yC  +  66  (2.10) 

where  a  is  the  coefficient  of  volume  expansion,  y  is  the  moisture  co¬ 
efficient  of  volume  expansion  and  6  is  the  coefficient  of  mechanical 
volume  expansion.  By  substitution  of  (2.10)  in  the  Doolittle  equation 
[28] 

log  a  =  -  2  303  f  '  7~  (2-11 ) 

oJ 

they  obtain  the  following  very  general  temperature,  moisture  and  pressure 
dependent  shift-factor 


log  a  = 


_ B  aAT  +  yC  +  66 

2.303  f  fQ  +  aAT  +  yC  +  60 


(2.12) 


Neither  equation  (2.8),  (2.9)  or  (2.12)  has  so  far  been  subject  to 
extensive  mechanical  testing  on  solid  polymers. 
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ii)  Single  Integral  Representations  Based  on  Thermodynamics 

In  order  not  to  obscure  the  main  issue  of  a  literature  review, 
we  shall  not  attempt  a  full  treatment  of  the  thermodynamic  foundations, 
but  rather  concentrate  on  the  hierarchy  of  different  continuum  thermo¬ 
dynamic  theories.  This  hierarchy,  shown  in  Fig.  2.3,  has  been  reviewed 
in  a  wel 1 -documented  survey  paper  by  K.  Hutter  [29].  Figure  2.3  shows 
that  the  main  concepts  which  can  be  utilized  to  develop  constitutive 
equations  for  inelastic  materials  are  the  rational  thermodynamics 
approach  and  the  state  variable  approach  (also  called  irreversible 
thermodynamics).  The  latter  will  be  discussed  in  detail  in  Chapter  3. 
Coleman  and  Noll  [30]  applied  functional  analysis  to  continuum  mechanics 
in  order  to  express  current  stress,  free  energy  and  entropy 
as  functionals  of  the  thermokinematic  history.  Their  analysis 
continues  to  influence  rational  thermodynamics  research  up  to  this  date. 

The  state  variable  approach  includes  certain  internal  variables 
in  order  to  represent  the  internal  state  of  the  material.  Constitutive 
equations  which  describe  the  evolution  of  the  internal  state  are  in¬ 
cluded  as  part  of  the  theory.  Onsager  [31]  published  the  concept  of 
internal  variables  in  thermodynamics.  This  formalism  was  used  by 
Biot  [32]  in  the  derivation  of  constitutive  equations  of  linear  visco¬ 
elasticity.  Ziegler  [33]  extended  this  work  to  include  nonlinear 
viscoelastic  materials  and  his  approach  has  been  used  by  Schapery  [34] 
and  Valanis  [35]  to  develop  specific  theories.  The  extension  to  in¬ 
clude  viscoplastic  materials  was  made  by  Perzyna  [36]. 
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Thermodynamics 


Ui 
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2.2.3  Mechanical  (Rheological)  Models 

Elastic  behavior  is  compared  to  that  of  a  massless  Hookean  spring 
and  the  viscous  behavior  to  that  of  a  Newtonian  dashpot.  The  success¬ 
ful  description  of  the  stress  strain  behavior  depends  on  the  proper 
choice  and  arrangement  of  the  elements  or  hypothetical  rheological 
units.  The  final  hardware-system  that  describes  the  behavior  of  a 
material  system  in  the  linear  range  can,  at  least  in  principle,  be 
extended  into  the  nonlinear  range  by  the  introduction  of  nonlinear 
springs  and  dashpots.  An  example  is  the  non! inearized  generalized 
Kelvin  element  whose  mathematical  equation  for  creep  can  be  written 
as , 

f9(o)  N  1 

e  =  f-  f,  (a)  +  — —  •  t  +  Mo)  l  p-  (1  -  exp  (-t/x. ))  (2.13) 

o  n  J  i=l  Li  1 

where  f. (a)  is  a  nonlinearizing  function  of  stress.  Bach  [37]  used 
this  form  in  his  studies  on  wood.  He  separated  the  total  creep  strain 
into  instantaneous,  delayed-el astic  and  flow  components.  The 
instantaneous  or  glassy  response  was  always  independent  of  stress 
(thus  f i (a )  =  1),  but  the  delayed  elastic  and  flow  compliances  were 
independent  of  stress  only  within  certain  limits  of  stress,  depending 
on  moisture  content  and  temperature. 

We  will  show  (Chapter  3)  that  Schapery's  thermodynamic  model,  at 
least  under  certain  circumstances,  can  be  related  to  a  spring-dashpot- 
hardware  system. 


14 


2.2.4  Power  Laws 


It  can  be  shown  that  the  use  of  discrete  spring-dashpot  models 


will  allow  us  to  obtain  an  accurate  representation  of  viscoelastic 
phenomena  over  a  large  time  scale  provided  we  use  a  sufficient  number 
of  elements  -typically  one  per  decade. 

The  idea  behind  the  power  law  representation  is  to  sacrifice 
accuracy  at  any  one  time  in  order  to  obtain  a  reasonable  representation 
over  the  complete  time  scale  between  glassy  (short  term)  behavior  and 
rubbery  (long  time)  behavior.  This  is  often  called  broad-band  approxi¬ 
mation.  Landel  and  Fedors  [38]  gave  a  form  for  spectral  repre¬ 
sentation  H(t)  guided  by  their  understanding  of  polymer  mechanics 

which  may  be  written  as, 
r  ^n 


H(t)  =  C 


T 


^  ) 


exp  (-  tq/t) 


(2.14) 


Clauser  and  Knauss  [39]  gave  examples  of  this  approximation  of  H(x) 
which  are  shown  in  Fig.  2.4.  These  results  lead  to  the  following 
broad-band  approximation  for  the  relaxation  modulus: 


Er(t)  =  E 


+ 

e 


(2.15) 


The  representation  is  seen  to  possess  four  arbitrary  constants, 
including  the  limiting  moduli  values.  The  exponent  n  gives  the  slope 
of  the  relaxation  curve  (Fig.  2.5)  through  the  transition  region 
between  glassy  and  rubbery  behavior  and  x  ,  for  isothermal  conditions, 
fixes  the  characteristic  relaxation  time. 


15 


Eq-Ee  tv  n 
H (r) :  — ^ — —  (-f  )  exp 
r/n)  '  T  ' 


" 1  7  -  V 


Recognizing  that  relaxation  and  retardation  are  basically 
reciprocal,  one  could  also  postulate  a  modified  power  law  for  the 
compliance  spectrum,  L(-t),  as 

f  1 n ' 

L(x )  =  C'  -T  exp(-x/T0)  (2.16) 

^  o' 

However,  the  integration  for  the  creep  compliance  does  not  give  a 
simple  form  and  it  is  more  direct  to  assume  a  reasonable  expression  as 


where  n  may  be  the  same  exDonent  as  in  the  relaxation  modulus  and  r' 

o 

is  chosen  for  the  best  fit  of  the  data.  A  special  case  results  when 
t  «  Tq  and  C  =  Dg  -  D  , 

n 

Mt>  -  D,  +  c  f|r|  (2.18) 

3  ^  V 

A  slightly  modified  form  can  be  written  as 

ft)n 

Dc(t)  -  Dg  +  C*  |yj  (2.19) 

where  =  A  U,  C*  =  C/An,  A  is  a  dimensionless  multiplication  factor, 
and  U  =  1  is  a  time  unit  selected  to  accommodate  data  reduction  (i.e. 

1  sec,  1  min,  1  hr,  etc.).  The  power  law  in  this  form  (Eq.  2.19)  is 
often  called  the  "Findley  power  law"  and  has  been  found  to  give  an 
accurate  representation  of  creep  behavior  of  an  extremely  broad  range 
of  amorphous,  crystaline,  crosslinked  glass  and  carbon-fiber  reinforced 
plastics  and  laminates.  The  equation  has  been  verified  mostly  in 
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tensile  tests,  but  in  a  few  cases  it  has  successfully  described 
behavior  in  compression  and  under  combined  tensile  and  shear  stresses 
as  wel 1 . 

Upon  conducting  a  screening  experiment  on  934  neat  resin  (at  room 
temperature),  we  identified  the  exponent,  n,  in  Eq.  2.19  as  the  most 
sensitive  indicator  of  experimental  errors.  In  addition,  n  was  shown 
to  be  dependent  on  the  duration  of  the  creep  test  as  shown  in  Fig.  2.6. 
Thus  any  short  time  creep  test  will  result  in  an  underestimation  of  n 
or  in  other  words  a  stable  evaluation  of  n  is  impossible  as  long  as  the 
flow  term  (C*tn)  is  too  small  relative  to  the  initial  compliance  D  . 

The  problem  is  easily  handled,  by  performing  a  small  additional  experi¬ 
ment,  which  consists  of  recording  the  recovery  strain  after  unloading. 
Creep  recovery  is  mathematically  given  by 

Sr  -  4S[(1  +  X)"  -  (X)"]  (2.20) 

t  -  1 1  n 

where  X  =  — ^ —  ,  as  =  a*C*t^  ,  and  t^  is  the  time  at  unloading. 

Eq.  2.20  can  be  very  easily  derived  using  the  power  law  for  creep  to¬ 
gether  with  the  superposition  principle.  Fig.  2.7  shows  that  the  n 
value  obtained  from  recovery  data  is  nearly  insensitive  to  the  duration 
of  the  recovery  experiment. 

The  creep  experiment  was  repeated  on  T300/934  carbon  epoxy  with 
the  load  perpendicular  to  the  fiber  direction.  The  result  plotted  in 
Fig.  2.6  is  compared  with  the  result  obtained  on  the  neat  934  resin 
(Fig.  2.8)  and  shows  that  the  composite  n-value  comes  within  5%  of 
the  resin  value  after  3000  min. 
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0.25 


Length  of  Creep  Test,  Log  t  (min) 


Fig.  2.6  The  Variation  of  Power  Law  Exponent  with 
Length  of  Creep  Test . 


Test  Duration  ,Log  t  (min) 


Fig.  2.7  The  variation  of  Power  Law  Exponent  with 
Length  of  Recovery  Test  . 
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Log  t  (min) 


Fig.  2.8  Comparison  of  Power  Law  Exponent 
Variation  with  Length  of  Creep  Test 


An  alternative  power-law  type  approach  was  developed  by  the 
Russian  school,  mainly  under  the  influence  of  the  ideas  published  by 
Yu.  N.  Rabotnov  [23].  His  formulation  was  based  on  the  viscoelastic 
deformation  of  polymeric  materials  which  can  be  described  by  an 
integral  equation  of  the  Bol tzmann-Vol terra  type 


e(t)  -  jr  [a ( t )  + 


,t 
J  o 


K(t  -  s)a(s)ds] 


(2.21) 


where  a(t)  and  e(t)  are  the  stresses  and  strains  in  a  uniaxial  stress 
field,  E  is  the  instantaneous  Young's  modulus,  and  K(t)  is  the  transient 
compliance  or  kernel  function.  Young's  modulus  can  be  determined  by 
the  basic  quasi  static  measurements,  while  the  kernel  must  include  the 
smallest  possible  number  of  parameters  that  need  to  be  determined 
experimentally.  To  satisfy  the  last  condition,  Rabotnov  approximated 
the  strain,  e(t),  with  a  power  law  and  proceeded  to  deduce  the  kernel 
function  K(t)  which  can  be  written  as. 


K(t)  =  Ata 


gn  tn(l+a) 
r'L(n+l  )(l+a)J 


(2.22) 


In  (2.22),  A,  3  and  a  are  material  parameters ,  r[(n+l ) (1+a)]  is  the 
gamma  function  of  the  argument  (n+l)(l+a),  where  1  +  a  >  0.  The  kernel 
parameters  can  be  determined  graphically  (when  a  full  deformation  curve 
is  obtained)  as  well  as  numerically  [19]. 

Functions  of  the  type  2.22  are  classified  as  exponential  functions 
of  fractional  order.  They  behave  as  an  Abel -type  singularity  for 
instantaneous  loading  and  like  an  exponential  for  the  succeeding  period 
of  time.  A  fairly  detailed  algebra  for  these  operators  has  been 
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developed  by  Rabotnov  and  his  collaborators. 

Equation  (2.21)  can  be  modified  for  the  case  of  creep  where  o(t) 
is  constant  and  becomes, 

e(t)  =  [i  +  f*  K(s )ds]  (2.23) 

-1  o 

For  small  t,  the  first  term  in  (2.22)  dominates,  i.e., 

“w-htkt  (2-24> 

By  substituting  (2.24)  into  (2.23),  we  obtain: 

6(t)  =  r [1  +  riky  tl+“]  (2-25) 

This  again  reduces  to  the  Findley  equation  (2.19) 

D(t)  -  DQ  +  C*tn 

where  C*  =  -(-['+n y  and  n  =  1  +  a.  It  is  clear  again  that  Eq.  2.19  is 
only  able  to  describe  a  limited  portion  of  the  full  deformation  curve. 
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Chapter  3 

SINGLE  INTEGRAL  REPRESENTATIONS  BASED 
ON  IRREVERSIBLE  THERMODYNAMICS 


3.1  Introduction 

A  single  integral  nonlinear  viscoelastic  constitutive  equation 
has  been  developed  by  Schapery  in  a  series  of  publications  [34,40-44]. 
This  chapter  reviews  his  thermodynamic  development  with  a  special  view 
to  developing  a  better  understanding  of  the  nature  of  the  four  non¬ 
linearizing  parameters  which  appear  as  well  as  their  physical  origin. 

The  name  irreversible  thermodynamics  is  well  established  but 
confusing,  since  it  suggests  geometric  irreversibil ity .  This  is  not 
generally  true  because  thermodynamic  irreversibility  relates  to  entropy 
production  while  a  deformation  accompanied  by  a  change  in  entropy  can 
be  geometrically  reversible  as  for  example  in  rubber  elasticity.  Thus 
a  name  like  "process  thermodynamics"  or  non-equilibrium  thermodynamics 
would  be  more  related  to  the  physics  of  a  problem  like  viscoelastic 
deformation. 

The  consistency  of  constitutive  equations  with  the  fundamental 
theorems  of  thermodynamics  has  been  discussed  by  Truesdell  and  Toupin 
[45].  No  restriction  on  constitutive  equations  is  given  by  the  first 
and  the  third  fundamental  theorems.  The  second  excludes  processes 
with  negative  entropy  production  and  thus  limits  the  field  of  constitu¬ 
tive  equations,  but  not  to  the  extent  desired.  It  was  demonstrated 
by  Biot  [46]  that  the  more  precise  statements  about  the  nature  of  the 
second  law,  which  were  made  by  Onsager  [47],  are  sufficient  to  treat 
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problems  in  linear  viscoelasticity. 


3.2  Gibbs  Free  Energy  Formulation  of  the  Schapery  Equations 

This  formulation  is  particularly  useful,  since  it  allows  the 
treatment  of  temperature  and  stress  tensor  components  as  independent 
variables.  The  Gibbs  free  energy  g  is  defined  as 

g  =  u  -  Ts  -  -  Qm  qm  (3.1 ) 

where  u  is  the  specific  internal  energy,  T  the  absolute  temperature, 
s  the  specific  entropy,  p  the  density,  q^  (m  =  l,2,...,k)  a  set  of  k 
state  variables  consisting  of  observed  variables  (strains),  and 
(m  =  l,2,...,k)  is  a  set  of  k  generalized  forces  defined  by  the  virtual 
work  condition  .  Eq.  (3.1)  shows  that  free  energy  can  be 

accumulated  by  increasing  internal  energy  by  decreasing  entropy  or 
through  a  potential  energy  loss  of  the  external  loading. 

A  mathematical  expression  for  the  strength  of  the  entropy  source 
as  a  result  of  "irreversibility"  has  to  be  derived  first.  We  assume 
that  the  classical  concept  of  entropy  is  extendable  into  non¬ 
equilibrium  situations.  This  leads  to  the  following  result  which  is 
called  the  fundamental  equation  of  the  system 

Sj  =  sI(u,qi )  (i  =  1 ,2,. . . ,n)  (3.2) 

The  index  I,  identifies  the  system  under  study.  The  state  variables 
q.j  (i  =  k+l,...,n)  represent  internal  degrees  of  freedom.*  These  will 
allow  us  to  describe  mi cro-structural  rearrangements  related  to 


*The  number  n  of  independent  thermodynamic  variables  depends  on 
the  nature  of  the  system. 
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straightening  and  relative  sTiding  of  long  chain  molecules.  These 
variables  are  thus  eventually  observable,  but  certainly  not  control¬ 
lable.  For  this  reason  (i  =  k+l,...,n)  =  0  since  any  force 
associated  with  an  internal  degree  of  freedom  would  provide  us  with  a 
means  to  control  this  internal  coordinate.  The  entropy-increment 
dSj  created  by  a  change  in  du  and  n  variations  of  dq^  is  measured  by 
the  total  derivative  of  the  fundamental  equation,  i.e. 


dsI  = 


du  +  l 

i=l 


-  dq . 

do .  , 


(3.3) 


where  q.1  means  all  coordinates  are  held  constant  except  q^ .  The  incre¬ 
ment  du  is  controlled  by  the  first  law  of  thermodynamics  through 


du  =  dh  +  l 
i  =  l 


Q,  dq. 


(3.4) 


Thus  according  to  (3.4)  any  increment  in  internal  energy  is  produced 
by  an  infinitesimal  amount  of  heat  (dh)  crossing  the  boundary  and/or 
an  infinitesimal  amount  of  mechanical  work  p-^  Q^.  dq^  done  on  the  system. 

The  incremental  process  is  reversible,  when  dq^  =  0.  Then, 
du  =  dh  and  according  to  the  second  law, 

dh  =  TdSj  (3.5) 

By  substituting  (3.5)  into  (3.3),  we  obtain 


In  other  words,  the  temperature  measures  the  sensitivity  of  the 
internal  energy  for  entropy  changes.  Thus  (3.3),  written  now  under 
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the  restrictions  of  a  reversible  process  becomes 


Eq.  (3.9)  is  a  result  known  in  classical  thermodynamics  as  Gibbs  equa¬ 
tion  and  describes  the  associated  reversible  process,  that  causes  the 
same  increment  of  internal  energy  du  as  the  irreversible  process 
described  by  Eq.  (3.4). 

By  elimination  of  du  between  (3.9)  and  (3.4),  we  obtain 


TdSj 


n 

dh  +  l 

i=l 


-1 

P 


(Q-j  -  Q-jR)  dq. 


(3.10) 


with  Q. (k+1 , . . . ,n)  =  0. 

Thus  the  entropy  of  a  mass  element  changes  with  time  for  two 
reasons.  First,  because  entropy  flows  into  the  volume  element  occupied 
by  this  mass  element,  second  because  there  is  an  entropy  source  due 
to  irreversible  phenomena  inside  the  volume  element.  This  entropy 
source  is  always  a  non-negative  quantity. 

Since  the  computational  procedure  needed  for  the  derivation  works 
only  on  isolated  systems,  we  now  create  such  a  system  by  the  artifice 
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of  adjoining  a  large  heat  reservoir  II  to  the  physical  system  I. 

Eq.  (3.11)  states  that  the  entropy  is  an  extensive  property,  i.e.  the 
entropy  of  the  system  is  the  sum  of  the  entropies  of  the  subsystems 

dSj+j j  =  dSj  +  dSjj  (3.11) 

and 

dsn  =  -  j1  (3.12) 

Eq.  (3.12)  states  that  the  reservoir  plays  the  role  of  a  reversible 
heat  source,  i.e.  heat  can  be  extracted  from  it  without  a  change  in  its 
temperature  T.  A  combination  of  (3.12),  (3.11),  and  (3.10)  leads  to 


'  5*  i=l  *< 


(3.13) 


D 

with  =  Q..  -  Q..  .  It  turns  out  that  the  strength  of  the  entropy 
source,  Sj+jj  has  a  very  simple  appearance.  It  is  a  sum  of  n  terms, 
each  being  a  product  of  a  flux  characterizing  an  irreversible  process 
and  a  quantity,  called  thermodynamic  force. 

In  order  to  proceed  from  here  we  need  to  calculate  the  rate  of 
approach  to  equilibrium.  This  is  only  possible  by  arguments  outside 
the  scope  of  continuum  mechanics.  The  entropy  source  is  explicitly  re¬ 
lated  to  the  various  irreversible  processes  that  occur  in  the  system 
through  a  linear  phenomenological  relation  between  forces  and  fluxes 
q.j  called  Onsager's  theorem. 


n 

=  l 

j=l 


b.  .(Q  ,T)o. 
ij  vym’  '  'j 


(3.14) 


This  fundamental  theorem  states  that  when  the  flow 


corresponding 
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to  the  irreversible  process  j,  is  influenced  by  the  force  X.,  cor¬ 
responding  to  the  irreversible  process  i,  then  flow  cm  is  then  also 
influenced  by  the  force  X.  through  the  same  influence  coefficients, 

J 

b.  .  (Qm,T) .  Eq.  (3.14)  can  be  rewritten  as 

Q.  =  Q.R  +  b. - (Q  ,T)q  .  (3 

ijwm’  'Mj  v 

It  is  shown  in  Appendix  B  that: 


3q. 


i  =  k+1 , . . .  ,n 


D 

The  generalized  force  Q..  associated  with  q^  through  the  Gibbs  equa¬ 
tion  (3.9),  is  derivable  from  the  Gibbs  free  energy  g  (Qm ,  qr,  T) 


qm  =  "  p  Iq^ 


m  =  1 ,2 , . . .  ,k 


The  problem  is  thus  reduced  to  the  solution  of  two  systems  of  equations. 


_  _  sg 

%,  -  -  0  aT 

m 


m  =  1 ,2 , . . .  ,k 


(3.16a) 


p  +  b  (Q  ,T)q 
3qr  rsv  m  /ys 


r  ,s  =  k+1  , . . .  ,n 


(3.16b) 


A  constitutive  assumption  on  the  dependence  of  the  Gibbs  free  energy 
on  Qm,  q^  and  T  has  to  be  made.  This  is  accomplished  by  a  Taylor  series 
expansion  around  the  equilibrium  state  such  that, 

pq  =  cn  +  c  Q  +  c  q  +id  Q  Q  +  d  Q  q 
y  R  m  ym  r  Mr  2  mn  ym  wn  rm  vm  Hr 

*?d,sVs  <3J7> 

where  c^,  c^  and  d^.  may  depend  on  temperature,  but  not  on  q^  and  Qm. 
Note  that  there  exists  a  distinct  relation  between  summation  index 
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and  physical  dimension  for  these  coefficients,  i.e. 


Coefficient 

Represents 

Summation  Range 

CR 

Energy 

cm 

Mechanical  strain 

m  =  1 ,2, . . .  ,k 

cr 

Thermodynamic  force 

r  =  k+1 , . . .  ,n 

dmn 

Mechanical  compliance 

n ,m  =  1 ,2, . . .  ,k 

Transforms  hidden  variables 
gr  into  mechanical  strain 
units 

r  =  k+1 , . . .  ,n ; 
m  =  1 ,2,. . . ,k 

drS 

Thermodynamic  modulus 

r,s  =  k+1 , . . .  ,n 

A  further  approximation  of  the  free  energy  is  found  by  an  expansion 
with  respect  to  the  uncontrollable  variables  only,  i.e.  we  can  be 
far  away  from  equilibrium  with  respect  to  Q  . 

09  "  pgR  +  cr  qr  +  I  drs  qr  qs  (3'18) 

g^,  and  d^$  can  be  functions  of  the  independent  variables  Qm  and  T. 
Eq.  (3.18)  must  reduce  to  (3.17)  under  sufficiently  small  forces, 
which  leads  to  the  relations 


Pdn  =  cn  +  c  Q  +id  Q  Q 
yR  R  m  wm  2  mn  wm  wn 


(3.19) 

(3.20) 

(3.21) 

The  various  terms  in  equation  (3.18)  have  the  following  physical 
significance; 


c 1  =  c  +  d  0 
r  r  rm  vm 


d'  =  d 
rs  rs 
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•  pg^  An  energy  term  that  acts  as  a  potential  for  the  time  inde¬ 

pendent  response.  The  dependence  of  cm  and  dmn  on  stress 
and  temperature  can  be  lumped  together  with  c^>  into  a 
single  nonlinearizing  parameter. 

•c'q  =(c  +d  Q)q  =  M  q  +  d  Q  q  This  term  is  de- 
r  r  r  rm  m  Mr  r  t  rm  m  t 

fined  by  the  energy  associated  with  a  simultaneous  action 

on  the  internal  coordinates  q  and  the  controllable 

r 

variables  Q  and  e. 
m 

•  d^$  q^  q$  This  term  is  defined  by  the  interaction  energy 

associated  with  the  simultaneous  action  on  internal  co¬ 
ordinates  q^  and  qs> 

The  above  items  which  make  up  equation  (3.18)  are  of  particular  interest 
as  they  define  the  basic  interaction  mechanisms  involved  in  three  dif¬ 
ferent  nonlinearizing  functions  which  appear  in  Schapery's  final  single 
integral  equation  as  developed  subsequently. 

Substitution  of  (3.18)  into  (3.16b)  gives 


c'  +  d1  (Q  ,T)q  +  b  (Q  ,T)  ^  =  0 
r  rs  m  Ms  rs  m  dt 


(3.22) 


In  order  to  reduce  (3.22)  to  a  system  of  n-k  differential  equations 
with  constant  coefficients,  the  further  assumption  is  made  that. 


brs«m-T>  “  VVT>brs 


in  which  brs  are  now  defined  as  viscosities.  This  equation  implies  that 
the  flow  is  Newtonian  but  with  stress  and  temperature  dependent 
viscosity.  The  implicit  stress  and  temperature  dependence  can  be  made 
explicit  through  a  shift  function  ^(Qp^T).  A  lot  of  experimental 
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evidence  exists  for  the  special  case  where  all  viscosities  have  a 


common  shift  factor  independent  of  stress.  In  such  cases  aD  becomes 
the  familiar  ay(T)  shift  factor  found  in  polymer  literature  [28]. 

The  following  form  for  d^$  is  guided  by  (3.23), 

drs(Qm'T)  =  aGWm'T)drs  r.s  -  k+1 ,n  (3.24) 

where  the  constant  symmetric  "thermodynamic"  modulus,  dJJs ,  is  now  de¬ 
fined  at  the  reference  temperature  TD  and  Q  =  0  as  indicated  schematical- 

k  m 

ly  in  Fig.  3.1  and  aG(Qm,T)  is  a  shift  function  as  defined  below. 

d«  =  d;s(°>TR)  =  drs<TR>  (3-25) 

It  follows  from  Eq.  (3.13),  (3.14)  and  (3.23)  that  a^  nonl inearizes 
the  entropy  production, 


SI+II  "  PT  aD^m’T)  bi j  qi  q; 


(3.26) 


On  the  other  hand  a^  nonl inearizes  the  third  term  in  the  Gibbs  free 
energy  expansion  (3.18)  and  therefore  pg  is  proportional  to  the  above 
quantities  as  indicated  by, 


pg  ~  2  W1)  drs  qr  qs  ^3'27'> 

We  cannot  make  an  a-priori  statement  on  the  dependence  of  ar  with 

b 

stress,  but  the  dependence  of  d'  (Q  )  at  a  certain  temperature  T  can 
increase,  decrease,  or  be  independent  of  the  force  Qm  (Fig.  3.2). 

A  constitutive  assumption  on  aG  has  been  developed  by  Schapery 
[48]  for  the  case  where  microcracking  and  flaw  growth  (in  the  opening 
mode)  are  the  dominant  softening  (damage)  mechanisms.  The  reciprocal, 
l/aG,  is  an  increasing  function  of  the  q-th  order  Lebesque  norm  of 
tensile  stress,  <}>',  as  specified  by. 
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=  1  + 


<f>' 


G(4>)  d<j> 


(3.28) 


where, 


rZ 


cr"  dC 


(3.29) 


G(<f>)  is  positive  and  reflects  the  statistical  distribution  of  damage. 

r 


The  quantity,  G ( )d<f> ,  can  be  found  from  mechanical  property  tests 


wherein  q  is  a  positive  constant. 

Extensions  of  Eq.  (3.29)  have  been  made  for  multiaxial  stress 
conditions  through  incorporation  of  the  first  and  second  invariants  of 
the  stress  tensor  (0  and  J^),  while  aging  and/or  rehealing  is  intro¬ 
duced  through  a  function  M(t)  such  that, 


1  a 

Mn  [A-j  e  +  A2/JJ  ]q  dc  q 


(3.30) 


Here,  A-j  and  A2  are  constants  which  satisfy  the  relation  A2  =  (1  - 

A-|)  and  when  they  are  both  positive  there  is  a  suppression  of  crack 
growth  by  (0  <  0). 

Due  to  the  simplifications  introduced  in  (3.23)  and  (3.24),  a 
system  of  n-k  differential  equations  with  constant  coefficients 
solvable  for  q^+p...,qn  now  remains.  These  together  with  the  k 
algebraic  equations  obtained  by  substituting  the  free  energy  expression 
(Eq.  3.18)  into  the  k  equations  of  (3.16a)  are, 

(n-k  equations)  (3.31 ) 


d9c 


aG  drs  «s  brs3T=-  S 


39 


qm  =  -  p 
m 


R 


3C' 


3Q„ 


> 

m 


(k-equations ) 


(3.32) 
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The  former  equations  (3.32)  can  be  rewritten  with  the  reduced 


time  4/  as, 

■R  h  d«s  _  -  cr 

^rs  qs  ^rs  di^  a^ 

where, 


(3.33) 


dijj  =  dt/a  and  a  =  an  I  a~ 

o  a  U  1  G 

The  final  system  of  equations  can  be  explicitly  written  as, 


S+1 

3Q1 

3ck+2 

'  aQ, 

K 

3Q1 
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(3.34) 


Remark  that  the  matrix  partition  on  the  upper  left  hand  corner  side 
has  a  maximum  dimension  of  6  x  (n-6)  while  the  two  lower  partitions 
are  square  matrices  with  dimensions  (n-6)  x  (n-6)  where  n  can  be  very 
large,  since  this  is  the  number  of  internal  coordinates,  needed  for  an 
adequate  description  of  the  inelastic  material  phenomena. 
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It  can  be  proven  that  a  simultaneous  diagonal ization  of  and 

rs 

£>rs  is  always  possible,  which  means  physically  that  a  transformed 
space  exists  where  the  n-k  internal  deformation  modes  do  not  interact. 
Thus , 


rs 


and 
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for 


(3.35) 
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d  q  +  b 

r  Mr  r  d^»  ; 

The  solution  of  (3.36)  is 


=  (1  -  e'+/TH 


(3.36) 


r  ' 

1 

c 1  ^ 
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dr 
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(3.37) 


with  =  b^/d  .  Eq.  (3.37)  has  a  fading-memory  Kelvin-type  kernel  and 
can  thus  be  seen  as  the  evolution  equation  for  a  Kelvin  unit  associated 
with  the  physical  process  represented  by  q^.  A  point  of  interest  here 
is  that  there  is  no  force  associated  with  the  (n-k)  functions  qr 
which  describe  the  evolution  of  the  internal  system.  It  appears  however 
that  these  internal  mechanisms  are  driven  by  forces  -  c^/a^  spiked 
into  the  Kelvin  unit  through  a  force  transfer  mechanism,  which  was 
mathematically  given  as  and  includes  a  thermodynamic  force  cr  due  to 
temperature  differences  (<j>  =  T  -  T^). 


c'  =  c  +  d  (Q  ,T)Q 
r  r  rnrxm’  'ym 


(3.20) 


The  transfer  function  (Eq.  3.20)  describes  how  forces  Qm  associated 
with  the  observed  coordinates  (q-,  ^ , . . .  .q^)  diffuse  down  into  the 
internal  material  structure  where  they  drive  the  hidden  coordinates 
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This  force  transfer  process  is  not  necessarily  linear  and  introduces 
a  third  nonlinear  function  Qm  through  the  equation 

c'  =  g  <f>  +  dR  Q  (Q  ,T)  (3.38) 

r  r  Y  rm  xnrMn 

where  <j>  is  T  -  TD  and  6  is  the  thermal  modulus. 
k  r 

Note  that  functions  d1  (which  show  up  in  the  third  term  of  the 

rs 

free  energy  expansion  (3.18))  and  drm  (which  show  up  in  the  second  term) 
represent  two  completely  different  nonlinear-physical  phenomena.  This 
reasoning  leads  to  the  two  distinct  functions  a^  and  Qm  respectively. 
Substitution  of  (3.37)  into  (3.31)  gives 


39c 


9C. 


°*m 


3Q 


m 


+  I 

r 


3Q, 


(1  -  e"^1^) 


m 


aP 
r  G 


(3.39) 


and  for  time  varying  forces. 
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By  substituting  (3.38)  in  (3.40), 
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np  Lv  np 

Aan  =  l  a<r)  (1  -  e“*/Tn 


(3.41) 


(3.42a) 

(3.42b) 
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(3.42c) 
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(3.42d) 


For  uniaxial  loading  where  q-j  is  the  tensile  strain,  e,  and  Q-j  is  the 
tensile  stress,  a,  Ea.  (3.41)  can  be  written  as 


r  -  39r  +  3Q]  ff*  -nr i  cl  d(Ql/aG}  h 
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(3.45a) 


91  =  3Qi 


(3.45b) 
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(3.45c) 


aa  ~  aD  1  aG 


(3.45d) 


By  way  of  summary,  the  nonlinear  parameters  gQ,  g-j ,  g2  and  a^ 
have  the  following  physical  significance. 

•  gQ  Non! inearizes  the  elastic  response  and  can  in  principle 

increase  or  decrease  with  stress.  The  former  is  described 
as  a  softening  phenomenon  while  the  latter  is  considered 


to  be  a  hardening  phenomenon. 
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g-|  Experimental  results  indicate  that  g-j  is  an  increasing 

function  of  stress  which  requires  Q-|  to  be  a  monotoni cal  ly 
increasing  function  of  Q-j  as  seen  in  Fig.  3.3.  This  means 
that  the  stresses  that  are  transferred  to  the  microstructure 
may  have  a  higher  intensity  than  the  external  force.  Note 
also  that  when  the  material  is  linear  viscoelastic, 
g-j  =  1  and  Q-j  is  a  linear  function  of  Q-j . 


Fig.  3.3  Dependence  of  Generalized  Force  Q-|  or  Q-| 


g2  This  factor  contains  two  nonlinearizing  functions,  Q-|  and 
a^.  It  can  be  seen  as  the  ratio  of  the  force  that  could 
be  spiked  into  Kelvin-type  microunits  taking  into  account 


that  the  moduli  of  the  springs  in  these  units  change  by  a 
factor  a^.  This  change  can  be  stiffening  (a^  >  1 )  or 
softeni ng  (ar  <  1  ) . 


a 

a 


Represents  a  stress  and  temperature  dependent  shift 
along  the  time  axis.  Here  ap(Qm’"0  decreases  with 
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temperature  and  with  increases  in  stress-induced  free 
volume.  Also,  a q ( Qm » T )  decreases  with  temperature 
but  an  increase  in  stress  would  produce  more  shift. 

The  Kelvin  type  kernel  aD  in  Eq.  (3.44)  can  be  approximated  by  a 
power  law  which  was  proved  to  be  valid  under  restricted  conditions  in 
Chapter  2, 

aD^-V)  =  l  D  (1  -  e'('M',)/Tn  (3.46) 

r 

or 

ADfW)  *  c(^*)n  (3.47) 

With  this  approximation  Eq.  (3.44)  becomes, 
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For  a  creep  test  a  =  aQ  H(t),  eq.  (3.48)  reduces  to, 


6  =  90  Do  °o  +  gl  g2  C 


J 

K  aJ 


(3.49) 


The  experimental  procedure  that  should  be  followed  to  obtain  the 
three  linear  parameters  (D  ,C,n)  and  the  four  nonlinearizing 
parameters  (gQ ,g-|  ,g£ ,a  )  has  been  described  in  detail  by  Lou  and 
Schapery  [49]  on  the  basis  of  a  graphical  shifting  procedure  based  on 
creep  and  creep-recovery  data.  A  computerized  procedure  that  avoids 
tedious  graphical  shifting  was  developed  during  the  current  research 
effort  by  Bertolotti  et  al .  [50]. 

The  general  equation  (3.44)  can  be  simplified  to  model  specific 
thermorheological  complex  materials.  For  sufficiently  low  stresses, 

-  Qm  =  Qm  =  a  which  means  that  g-|  =  1  and  ~  "*/aQ-  Also,  a^  and  aQ 
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are  only  functions  of  temperature  Up/a^  =  ay) . 
to 
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dx 
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Eq.  (3.44)  now  reduces 

(3.50) 


where  log  ip  =  log  t  -  log  a-p.  For  a  -  aQ  H(t)  this  becomes 


Dt  =  Dj(T)  +  ADU)/aG(T)  (3.51) 

with  Dj  =  e  |  aQ  and  D j (T)  =  gQ(T)  DQ.  Thus, 

log  (Dj  -  Dj)  =  log  aD  -  log  a^  (3.52) 


Eq.  (3.52)  shows  that  master  curves  can  be  obtained  by  rigid  horizontal 

1 1  og  a-j- 1  and  vertical  [  1  og  a^ |  shifts. 

Eq.  (3.52)  suggest  that  the  machinery  needed  for  accelerated 

characterization  can  only  be  effective  when  the  vertical  shift  ag(T) 

can  be  separated  from  the  horizontal  shift  a-|-(T).  This  is  only 

possible  with  additional  transient  temperature  tests  as  suggested 

from  Eq.  (3.41).  Some  actual  data  were  obtained  by  Schapery  and  Martin 

[51]  and  very  recently  by  Weitsman  [52].  In  the  former,  a  theoretical 

form  for  ar(T)  which  applies  below  and  in  the  neighborhood  of  the 
b 

glass  transition  temperature  Tg  was  given  as, 


a  =  y  W  (1  +  z) _  (3.53) 

G  exp  [y2/(l  +  z)]  -  1 

where  Y  =  aH/RTr  and  z  =  (T  -  TR)/TR.  Eq.  (3.53)  is  based  on  kinetic 
theory  wherein  aH  is  the  constant  activation  enthalpy  for  hole  forma¬ 
tion.  Simple  forms  for  a^  on  the  basis  of  rubber  elasticity  are  well 
established  at  and  above  the  Tg. 
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The 

material s 


results  discussed  herein  have  been  applied  to  composite 
by  Schapery  and  his  co-workers  [49,53]. 
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Chapter  4 
DELAYED  FAILURE 


4.1  Introducti on 

Failure  is  most  often  treated  as  a  separate  issue  from  the  deter¬ 
mination  of  modulus  properties  of  materials.  In  fact,  most  failure  laws 
are  derived  empirically  from  observations  related  to  a  catastrophic 
event  such  as  yielding  or  rupture.  As  a  result,  a  great  deal  of 
testing  and  data  analysis  is  necessary  to  establish  an  appropriate 
failure  law  for  a  material.  On  the  other  hand,  modulus  or  constitutive 
laws  are  derived  by  more  rational  means  of  relating  deformations  to  the 
forces  which  produce  them.  For  this  reason,  often  much  less  testing  is 
necessary  to  define  a  constitutive  law  for  a  material  especially  if 
deformations  do  not  depart  from  the  elastic  or  reversable  deformation 
range  for  a  material . 

Failure,  however  it  is  defined,  should  be  a  part  of  a  complete 
constitutive  description  of  a  material.  In  other  words,  the  key  to 
dealing  effectively  with  failure  lies  in  treating  its  behavior  as  a 
termination  of  a  nonlinear  viscoelastic  process.  Perhaps,  for  this 
reason,  a  number  of  investigators  have  suggested  that  modulus  and 
strength  laws  should  be  related  to  each  other  for  viscoelastic  materials 
[4,33,54]. 

The  concept  of  distortional  energy  as  a  measure  for  the  criticality 
of  a  given  state  of  stress  in  an  elastic  material  (von  Mises  criterion) 
cannot  be  carried  over  directly  as  a  governing  criterion  for  failure 
in  viscoelastic  materials,  mainly  because  viscoelastic  deformation 
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involves  dissipative  mechanisms.  Thus  at  any  time  the  energy  balance 
can  be  written  as: 

Total  deformation  =  Free  (stored)  +  Dissipated 
energy  energy  energy 

In  case  of  a  perfectly  elastic/viscoplastic  material,  inelastic 
deformations  occur  only  after  a  limiting  value  of  stress  has  been 
exceeded.  Owing  to  this  assumption,  such  a  yield  criterion  does  not 
differ  from  the  known  criteria  of  the  inviscid  plasticity  theory.  The 
most  appropriate  form  on  the  latter  is  often  based  on  the  energy  approach 
(Huber-Henkey-von  Mises)  in  which  a  critical  value  of  the  conserved 
(accumulated)  distortion  energy  defines  the  transition  from  the  elastic 
to  the  inelastic  state. 

In  case  of  an  elastic-viscoplastic  material,  the  rheological 
response  occurs  from  the  beginning  of  the  loading  process.  Under  such 
conditions  it  seems  appropriate  to  take  into  account  the  accumulated 
energy  W^,  as  well  as  the  dissipated  power  Wp  where  both  quantities  are 
regarded  to  be  responsible  for  the  behavior  of  the  viscoelastic  material. 
The  yield  condition  is  then  given  by  a  function  with  two  arguments 
such  that 

f(wE,wD)  =  0 

where  the  specific  form  is  based  on  experimental  evidence.  The  simplest 
possible  form  is  the  linear  combination 

«E  WE  +  «D  CliD  =  k2 

with 

0  1  {^e,CD}  -  1 
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With  Cq  =  0  we  assume  that  the  dissipated  energy  has  no  effect  on 
the  yield  criterion.  This  corresponds  to  Reiner  and  Wei ssenberg ' s 
point  of  view.  Reiner  and  Weissenberg  [55]  suggest  that  the  energy 
storage  capacity  is  responsible  for  the  transition  from  viscoelastic 
response  to  yield  in  ductile  materials  or  fracture  in  brittle  ones. 

They  assume  that  a  threshold  value  of  the  distortional  free  energy, 
called  the  resilience  of  the  material,  is  the  governing  quantity.  It 
is  clear  that  the  Reiner-Wei ssenberg  criterion  applied  to  a  material 
with  zero  dissipation  (elastic  material)  becomes  identical  to  the  von 
Mises  criterion.  When  applied  to  a  viscoelastic  material,  however,  we 
have  to  keep  track  of  the  free  energy  which  is  history  dependent.  If 
the  mechanisms  through  which  total  deformation  energy  is  transformed 
into  dissipated  energy  are  so  activated  that  no  free  energy  can  be 
accumulated,  there  is  practically  no  limit  to  the  amount  of  deformation 
energy  which  can  be  applied  without  impairing  the  strength.  That  is, 
forces  up  to  a  certain  magnitude  can  be  applied  for  any  length  of  time 
without  leading  to  rupture.  If  it  turns  out  that  the  material  cannot 
accommodate  this  energy  redistribution  fast  enough  then  the  material 
will  store  energy  initially  but  after  a  period  of  time  failure  will 
occur.  The  failure  process  is  therefore  delayed  which  limits  the  life 
of  a  given  structure  to  a  finite  value.  The  instant  of  yielding  is 
thus  clearly  dependent  on  the  final  outcome  of  a  conjection  between 
deviatoric  free  energy  and  dissipated  energy.  Thus  the  effect  of  the 
strain  history  on  the  delayed  yielding  phenomenon  follows  from  this 
model  in  a  natural  way.  In  this  respect,  it  should  be  noted  that  both 
Naghdi  and  Murch  [56]  and  Crochet  [57]  have  used  a  modified  von  Mises 
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criteria  for  viscoelastic  materials  by  assuming,  ab  initio,  that  the 
radius  of  the  yield  surface  Y  does  depend  upon  the  strain  history, 
through  a  Euclidian  metric,  x>  in  strain  space.  Their  formulation 
can  be  expressed  as , 


ij  sij  -  y2  -  0  •  Y  ■  v<*> 

(4.1) 

V.  .  eE.)(6V.  -  6E.)]'/2 

1J  IJ  1J  IJ 

(4.2) 

where  x  represents  the  distance  between  the  purely  elastic  strain 
response  e^.  and  the  strain  state  corresponding  to  yield  eY.  as  is 

1  J  *  vJ 

geometrically  shown  in  Fig.  4.1.  The  functional  relationship  between 


Fig.  4.1  Geometric  interpretation  of  the  Euclidian  metric  x  (as  de¬ 
fined  by  Crochet). 

Y  and  x  for  one-dimensional  behavior  was  given  by  Crochet  [57]  to  be, 

Y(t)  =  A  +  Be_Cx  (4.3) 

where  A,  B  and  C  are  material  parameters.  The  radius  of  the  yield  sur¬ 
face  is  thus  a  continuously  decreasing  function  of  the  "distance"  between 
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two  states  of  strain  (e^  •  and  ev  . ) . 

An  extension  of  the  Naghdi  and  Murch-Crochet  approach  to  fiber 
reinforced  material  systems  is  possible  when  we  can  relate  the  stress 
in  the  matrix  phase  {a>m  to  the  applied  external  stress  {a}.  The 
introduction  of  stress  concentration  factors,  which  relate  both  these 
stress-fields  is  suggested  by  Hill  [58] 


{a}  =  B  {a} 

m  m 


(4.4) 


and  equation  (4.1)  becomes 


f  =  {a}T[Bm]{a}  -  Y2  =  0 


(4.5) 


in  which  Y  is  the  tension  yield  stress  in  the  matrix,  as  defined  by 
eq.  (4.3).  It  should  be  noted,  however,  that  eq.  (4.5)  only  describes 
initial  yielding.  The  incorporation  of  a  hardening  rule,  together  with 
a  flow  rule  should  be  considered  in  order  to  model  accurately  the 
phenomena  after  initial  yielding  up  to  final  fracture. 

The  extension  of  the  Reiner-Wei ssenberg  approach  to  a  general 
anisotropic  material  system  can  be  done  quite  easily  when  we  use  the 
general  definition  of  distortional  energy,  as  defined  by  von  Mises  [59] 

Hdist  5  7  <sn  4  +  s22  4  +  +  s66  Tz>  +  S1 2  Vy  +  "■ 


+  s14  Vx  +  •"  +  s66  Vz 

_  1  ^SH  +  S1 2  +  S1 3^  +  tz^s61  +  $62  +  s63^' 
2  sn  +  s22  +  s33  +  2^S12  +  s23  +  S31  ^ 


(4.6) 


The  advantage  of  this  approach  is  that  the  onset  of  yielding  or  failure 

is  controlled  by  only  one  material  constant  W, .  , ,  while  the  former 

di  st 
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approach  (Naghdi  et  al  .<■)  introduces  at  least  three  parameters.  It  has 
been  argued  by  Pagano  [60],  on  the  other  hand,  that  the  distortional 
energy,  as  defined  by  von  Mises,  has  only  physical  significance  in  the 
special  case  of  an  isotropic  material-  The  full  implications  of 
Pagano's  comments  are  unknown  at  the  present  time.  A  possible  alter¬ 
native  way  to  extend  the  Reiner-Weissenberg  criterion,  would  be  to 
evaluate  the  free  energy  in  the  matrix  phase.  This  could  also  be  done 
through  Hill's  concentration  factors. 

The  following  is  a  brief  review  of  the  Reiner-Weissenberg  criterion 
which  follows  Briiller's  [61-68]  extensive  investigation  along  these 
1 i nes . 

4.1.1  Free  Energy  Accumulation  for  a  Linear  Kelvin  Chain  in  Creep 

For  a  single  Kelvin  unit  in  series  with  a  spring  as  shown  in 
Fig.  4.2, 


/////// 


€,  (t) 


Fig.  4.2  Linear  Kelvin  model 
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where 


e  =  e  +  e, 

g  i 


e  =  D  a 

g  go 


e1  =  D1  ao^  "  exp  ("t/xl  )  ] 


(4.7) 


(4.8) 


e(t)  =  [D  +  D^l  -  exp  (-t/^JJc 


(4.9) 


The  stress,  a^,  in  the  spring  with  compliance  D-j  is, 


as  =  D7  =  ao  [1  "  exp  (-t/T!  ^ 


(4.10) 


The  stress,  a^,  in  the  dashpot  is  then, 
ad  =  ao  exp  ^"t/Tl) 


(4.11) 


A  physical  i nterpretation  of  free  energy  would  be  the  energy  stored  in 
the  springs  and  D-j , 

Springs  =  °o  2^  +  !  °s  ®s  dt 

J0 


9  D  9  D-j  ✓ 

=  0Q  2^-  +  aQ  2~  [1  -  exp  (-t/T-j)]' 


(4.12) 


In  the  same  way  we  can  calculate  the  energy  dissipation  in  the  dashpot, 


dashpot 


CTd  €d  dt 


[1  -  exp  (-2t/x1 )] 


(4.13) 


The  total  energy  in  the  three  parameter  system  is  given  by, 
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Jtotal  ^springs  +  ^dashpot 

9  FDn 

=  oQ  ^  +  Di  -  exP 


!-t/x1 ) 


(4.14) 


The  extension  of  Eq.  (4.12),  (4.13)  and  (4.14)  to  a  n-element  Kelvin 
unit  with  a  free  spring  is  straightforward. 
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9  [On  N  Di 

=  a2  p5-  +  l  J-  (1  -  exp  (-t/T.))‘ 

0  2  ih  2  1 


(4.15) 


N  D. 


dashpots  o  >-| 


l  21  (1  -  exp  (-2t/x.)) 


(4.16) 


“total  =  °o  ^  j,  Di  0  '  6XP  ('t/T1>> 


(4.17) 


By  taking  the  limits  for  t  ->  »  of  these  equations,  we  obtain 


1  i  m  W_ 


spri  ngs 


,  TD  N  D.l 
2  _3.+  y  _L 

°  L2  i=l  2 


(4.18) 


2  N  D. 
^dashpots  "  °o  ^  2 


(4.19) 
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Thus  half  the  work  done  by  the  external  forces  goes  to  increase  the 
free  energy  while  the  other  half  is  dissipated. 


4.1.2  Power  Law  Approximations  for  Free  Energy 

The  free  energy  expression  in  the  form  of  Eq.  (4.15)  forces  us  to 


identify  the  2N  material  constants  D.  and  x . .  This  can  be  avoided 
J  1  1 


Note  on  the  Zhurkov  Criterion 

The  existence  of  a  strength  limit  as  a  real  physical  characteris¬ 
tic  of  a  material  was  seriously  doubted  by  Zhurkov  [69]  due  to  experi¬ 
mental  evidence  that  solids  may  also  fracture  at  stresses  below  the 
strength  limit.  He  set  up  systematic  and  careful  measurements  of  the 
lifetime  of  solids  under  a  constant  tensile  stress.  Both  tensile 
stress  and  temperature  were  widely  varied  among  experiments.  More  than 


a  hundred  substances  of"all  principal  types  of  solids  were  investigated — 
metals,  alloys,  glasses,  polymers,  crystals.  It  turned  out  that  the 
time  to  fracture  possessed  a  surprisingly  uniform  dependence  when 
suitably  expressed  in  terms  of  the  stress  (a)  and  temperature  (T). 
Analytical  treatment  of  this  vast  amount  of  data  led  him  to  the  uni¬ 
versal  formula  for  the  temperature  and  stress  dependence  of  lifetime, 
x,  (except  for  very  small  stresses), 

x  —  XQ  exp  [(U0  -  y a  )/kT]  (4.26) 

where  k  is  the  Boltzmann  constant,  tq  is  a  constant  on  the  order  of 

1  3 

the  molecular  oscillation  period  of  10  sec.,  UQ  is  a  constant  for 
each  substance  regardless  of  its  structure  and  treatment,  and  y  depends 
on  the  previous  treatment  of  the  substance  and  varies  over  a  wide 
range  for  different  materials. 

Equation  (4.26)  can  be  formally  incorporated  into  a  continuum 
mechanics  analysis  in  a  modified  form  to  account  for  time  varying 
stresses  and/or  temperatures. 


Roylance  and  Wang  [70]  used  Eq.  (4.27)  in  a  finite  element  code  where 
the  current  value  of  the  above  integral  was  computed  at  each  node.  The 
time  and  location  of  fracture  was  determined  when  the  integral  value 
reached  unity  at  any  node. 

Williams  and  his  collaborators  [71,72]  disagree  with  Eq.  (4.26). 
They  argue  that  the  rupture  process  is  much  more  detailed  than  re¬ 
flected  by  this  equation.  Their  experiments  on  stress  history  induced 
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rupture  led  them  to  the  conclusion  that  "Zhurkov's  impressive  findings 
appear  to  be  the  somewhat  fortuitous  result  of  a  particular  test 
method. " 

When  Eq.  (4.27)  is  applied  to  such  highly  heterogeneous  material 
systems  as  composites,  the  question  of  whether  or  net  fracture  can  be 
regarded  as  a  thermally  activated  process  of  damage  accumulation  is  far 
from  settled.  It  is  also  obvious  that  x  ,  U  and  y  lose  their  physical 
meaning.  Regel  et  al .  [73]  suggest  a  life  of  mixture  type  predictions 
for  these  variables. 
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Chapter  5 

RESIN  CHARACTERIZATION:  EXPERIMENTAL  RESULTS 
AND  DISCUSSION 


5.1  Introduction 

In  order  to  understand  experimentally  observed  discolorations  of 
the  Fiberite  934  neat  resin  due  to  postcuring  and  some  anomalities  of 
the  creep  response  in  the  temperature  region,  200-250°F  (93.3-121 °C) , 
it  was  decided  to  study  details  of  the  relationship  between  resin 
chemistry,  various  postcuring  procedures,  and  mechanical  properties. 
General  resin  information  is  summarized  in  5.2  while  experiments  and 
analysis  related  to  the  mechanical  characterization  of  the  resin  are 
reported  in  5.3. 

5.2  General  Resin  Information 

5.2.1  Resin  Chemistry  and  Curing  Procedure 

Many  composites,  especially  carbon-fiber  based  composites,  utilize 
an  epoxy  resin  matrix.  From  the  hundreds  of  different  epoxies  that 
are  available,  only  a  small  number  {-  10)  are  used  for  composites.  The 
epoxy  of  interest  in  this  study  is  a  polyfunctional  epoxide  based  on 
the  tetraglycidyl  derivative  of  methyl enedi anil ine  (TGDDM),  available 
from  Ciba-Geigy  under  the  tradename  MY  720.  This  epoxide,  which  is  now 
the  basis  of  virtually  all  aerospace  matrices  for  carbon  fibre  both 
here  and  abroad,  is  shown  in  Fig.  5.1a.  The  Ciba-Geigy  Eporal  hardener 
used  in  conjunction  with  MY  720  is  4,4'  di ami  nodi  phenyl  sulfane  (DDS) 
Ciba-Geigy  Eporal  and  is  shown  in  Fig.  5.1b. 
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a)  Tetraglyci dyl  4,4'  di ami  nodi  phenyl  methane  epoxy  TGDDM 
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b)  4,4'  di ami  nodi  phenyl  sulfane  (DDS) 

Fig.  5.1  Chemical  Structure  for  Epoxy  Resins 

Because  the  DDS  cure  is  sluggish,  even  at  moderately  high  tempera¬ 
tures,  a  boron  trifluoride  catalyst  is  used  ( B )  in  order  to  produce 
a  more  manageable  cure  cycle.  Despite  the  fact  that  these  ionic 
curing  agents  often  have  poor  humidity  resistance,  they  are  used  as  a 
compromise  between  mechanical  properties  and  manufacturing  convenience. 
The  final  epoxy  network  structure  depends  not  only  on  the  chemistry  of 
the  system,  but  also  on  the  initial  epoxy  versus  curing  agent  ratio 
and  cure  conditions.  At  the  beginning  of  this  research  program,  we  ob¬ 
tained  three  different  panels  manufactured  with  three  different  curing 
procedures.  Cure  details  are  given  in  Appendix  A. 
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It  is  usually  assumed  that  the  properties  of  the  cured  resin  in 
a  composite  are  the  same  as  those  of  the  bulk  material.  However, 
McCullough  [74],  Hancox  [75]  and  Yang,  Carlsson  and  Sternstein  [76] 
suggest  otherwise.  The  reasons  mentioned  by  these  authors  are: 


•  The  structure  in  a  small  volume  of  resin  may  be  more  ordered 
than  in  the  bulk  material,  because  the  exotherm  in  the 
composite  will  be  limited  to  the  overall  curing  temperature 
by  the  presence  of  good  conducting  fibers. 

•  The  small  distance  between  the  fibers  may  provide  steric 
hindrance  and  inhibit  epoxy-amine  cross-link  reactions 
during  the  latter  stages  of  cure  and  limit  the  overall 
achievable  cross-link  density. 

•  Due  to  the  large  amount  of  surface  area  of  the  fibers, 
adsorption  of  epoxy  segments  at  the  fiber  surface  would 
cause  a  tighter  network. 

•  A  state  of  negative  hydrostatic  stress  could  be  imposed  on 
the  epoxy  due  to  the  fibers,  restricting  the  volume  expan¬ 
sion  of  epoxy  during  heating. 


No  convincing  evidence  has  been  produced  yet  (to  this  author's  knowledge) 
that  one  or  a  combination  of  these  mechanisms  leads  to  the  resin  in  the 
composite  that  is  thermorheological ly  different  from  the  neat  resin. 

As  a  polymer  is  cooled  at  a  finite  rate  through  the  glass  transi¬ 
tion  into  the  glassy  state,  a  sudden  increase  in  modulus  parallels-  a 
rapid  decrease  in  molecular  mobility.  The  process  of  molecular  re¬ 
arrangements  operates  on  its  own  internal  time  scale.  Equilibrium  can 
only  be  achieved  if  the  cooling  rate  is  so  slow  that  the  molecular  re¬ 
arrangements  are  given  sufficient  time  to  operate.  If  the  cooling  rate 
is  too  fast,  these  rearrangements  are  essentially  quenched  into  a  non¬ 
equilibrium  state,  i.e.,  a  state  characterized  by  excess  entropy, 
enthalpy  and  volume.  This  gives  rise  to  a  driving  force  towards 
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equilibrium,  i.e.,  the  system  densities  till  the  free  volume,  excess 
entropy  and  excess  enthalpy  all  reach  zero.  This  process  can  operate 
on  a  scale  of  days,  months  or  years  and  is  known  as  aging. 

Experimental  equilibrium  has  been  unexpectedly  difficult  to 
achieve.  Spontaneous  changes  in  the  properties  of  glassy  polymers  have 
been  observed  long  after  vitrification. 

Excess  entropy  cannot  be  measured  directly  but  excess  enthalpy 
is  easily  measured  by  calorimetric  methods  (DSC).  The  latter  show  a 
decrease  in  enthalpy  with  time  in  a  similar  fashion  as  the  time  dependent 
mechanical  properties  as  reported  by  Petrie  [77,78].  As  indicated  in 
Fig.  5.2,  she  did  not  find  a  one  to  one  correspondence  to  specific 
volume  behavior. 

In  our  experiments,  we  tried  to  achieve  this  equilibrium  state  by 
a  postcuring  of  350°F  ±  10°F  (176°C±5°C)  for  4  hr  ±  15  min.  followed 
by  a  slow,  controlled  cooldown  at  a  rate  of  5°F/hr  (2 ,5°C/hr)(see  Fig. 
5.3). 

5.2.2  Thermal  Transitions  in  the  Epoxy  System 

Two  distinct  transition  temperatures  at  140°F  (60°C)  and  356°F 
(180°C)  respectively  have  been  identified  by  Yeow  and  Brinson  [4]  dur¬ 
ing  their  experiments  on  T300/934  graphite  epoxy.  The  latter  was  the 
primary  glass-transition  temperature,  T  ,  while  the  former  was  a 
secondary  transition.  Their  measurements  were  based  on  detecting  dis¬ 
continuities  in  the  coefficient  of  thermal  expansion  which  were  measured 
using  electrical  strain  gages. 
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Fig.  5.2  Shear  Storage  Modulus,  6,',  Logarithmic  Decrement,  a,  and 
Enthalpy  Relaxation,  aH,  of  Quenched  Atactic  PS  vs. 
Annealing  Time  at  92°C.  G-j  and  A  Were  Measured  at  0.6  Hz. 
Petrie  [77,78]. 


Time  (  hr ) 


Fig.  5.3  Postcure  Cycle  Used  for  934-Resin  and  for  T300/934  Composite. 
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Temperature  (  deg  C  ) 


Keenan,  Seferis  and  Quilivan  [79]  used  oscillating  mechanical 
fields  to  characterize  the  dynamic  mechanical  behavior  of  an  epoxy 
system.  Their  results  are  shown  in  Fig.  5.4  with  an  w-transi tion 
centered  around  212°F  (100°C)  and  an  a-transition  at  464°F  (240°C). 
However,  these  definitions  of  Tg  suffer  from  the  fact  that  they  are 
based  on  dynamic  mechanical  quantities  which  are  quite  sensitive  to  the 
frequency  of  experimentation  as  well  as  the  sample  heating  rate.*  Con¬ 
sequently,  the  Tg  values  defined  from  dynamic  mechanical  experiments 
are  invariably  higher  than  those  obtained  by  a  more  conventional 
method,  i.e.,  dilatomeric  experiments. 


Temperature  (°C) 

a)  Storage  Modulus 


Temperature  (°C) 

b)  Loss  Modulus 


Fig.  5.4  Storage  Modulus  and  Loss  Modulus  after  Keenan  et  al .  [79]. 


*The  conclusions  obtained  by  Keenan  et  al .  can  be  carried  over 
to  the  934  resin,  since  this  particular  member  of  the  TGDDM-DDS 
epoxy  family  has  the  same  basic  network.  Any  difference  should  be 
one  of  degree,  not  in  kind. 
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5.2.3  Environmental  Stability 

1 )  Thermal  Stability 

Different  IR  analysis  of  an  epoxy  postcured  above  350°F  (180°C) 
reveals  an  additional  broad  sorption  band  in  the  5.5  -  14.3  y  region 
according  to  Morgan  [80].  The  intensity  of  this  sorption  increased  with 
exposure  time  and  increasing  temperature.  A  strong  IR  sorption  at 
5.85  y  appears  during  the  initial  stages  of  the  degradation  process, 
which  indicates  the  formation  of  carbonyl  groups.  Such  groups  should 
result  from  oxidation  of  unreacted  epoxide  rings  to  form  a-hydroxy 
aldehyde  and  carboxylic  acid  [81]. 

Degradation  causes  a  discoloration  of  the  934  epoxy  as  indicated 
in  Fig.  5.5  and  may  indicate  a  significant  change  of  behavior  in  the 
glassy  region  due  to  postcuring.  We  observed  that  heating  above  356°F 
(180°C)  for  36  hr  seemed  to  cause  embrittlement,  i.e.  decreased  tensile 
strength  and  ultimate  elongation  accompanied  by  a  slightly  increased 
modulus  (Fig.  5.6).  In  addition,  an  increased  moisture  absorption 
rate  (Fig.  5.7)  was  found.  This  is  most  likely  due  to  a  thermal  oxi¬ 
dation  phenomenon.  Figure  5.8  shows  the  storage  modulus  E'  and  the 
damping  ratio  6  for  the  same  resin. 

2)  Moisture  Absorption 

One  consequence  of  diffused  water  is  that  it  lowers  the  mechanical 
rigidity  of  the  resin.  In  extreme  cases  the  effect  is  to  lower  the 
glass  transition  temperature  by  as  much  as  212°F  (100°C).  The  effect 
of  moisture  on  the  glass  transition  temperature  and  the  rates  of 
absorption  and  desorption  for  different  resins  (including  the  Fiberite 


59 


As  delivered 


4  hours  at  350°  F  (  I77°C), 
slow  cooled  ,  5  °  F 
(  2. 5  °  C  )  per  hour . 


15  minutes  at  approximately 
10-15°  F  (  5- 8°  C  )  above  Tg  , 
then  quenched  down  to  room 
temperature . 


Short  overheating  at  500  °F 
(260°  C  ) ,  then  quenched  down 
to  room  temperature. 


36  hours  at  380°  F  ( 193° C) , 
slow  cooled  ,  5  °  F 
( 2.5  °  C  )  per  hour . 


Fig  5.5  Effect  of  Cure  Cycle  on  Resin  Color 
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934)  has  been  studied  by  Deisai  and  Whiteside  [82].  Figure  5.9  shows 
the  relationship  between  equilibrium  moisture  concentration  and  primary 
glass  transition  temperature  for  several  epoxy  resins.  Figure  5.10 
shows  the  effect  of  the  humidity  in  the  test  environment  on  the 
equilibrium  moisture  content  in  the  same  resins.  Figure  5.10  shows 
that  a  prolonged  postcure  period  could  adversely  affect  the  moisture 
absorption  rates,  indicating  a  damaged  matrix.  The  specimens  were 
allowed  to  pick  up  moisture  in  the  lab  [75-80°F  (24-27°C)  and  ~  45%  RH]. 
Figure  5.11  shows  the  influence  of  moisture  on  the  w-transi tion .  It 
is  seen  that  moisture  has  a  strong  effect  especially  on  the  width  of  the 
(jj-transition. 


Moisture  Content,  %(wt) 


Fig.  5.9  Effect  of  Moisture  Content  on  [82]. 
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The  specimens  we  used  for  creep  and  creep-recovery  lost  approxi¬ 
mately  1.5%  of  their  total  weight  during  postcuring.  For  this  reason, 
the  following  precautions  were  taken  to  control  moisture  absorption: 

•  The  specimens  were  dessicated  immediately  after  postcuring. 

(The  relative  humidity  in  the  dessicator  was  monitored  and 
kept  below  20%. ) 

•  All  specimens  were  tested  within  a  month  after  postcuring  when 
possible. 

•  The  specimens  tested  at  room  temperature  were  surrounded  with 
dessicant  that  was  kept  in  the  oven  space. 

In  conclusion,  our  own  moisture  absorption  experiments,  together  with 

a  careful  study  of  the  literature  available,  let  us  believe  that  the 

precautions  taken  were  sufficient  to  keep  the  moisture  level  well  below 

.3%  WT,  a  level  that  is  believed  to  be  negligible. 

5 . 3  Results  and  Discussion  on  934  Resin 

Two  934-resin  panels  were  fabricated  by  Graftech  as  detailed  in 
Appendix  A.  These  panels,  labeled  I  and  II,  were  rectangular  plates 
of  approximately  25"  x  20"  x  .12",  while  one  thick  plate  of  10"  x  10"  x 
.39"  was  made  by  Fiberite.  The  panels  were  cut  into  8"  x  1" 
rectangular  specimens  and  labeled  according  to  batch  and  position.  In 
order  to  have  as  homogeneous  a  set  of  specimens  as  possible,  each  was 
subjected  to  a  postcure  cycle  of  350°F  (176°C)  for  4  h  as  discussed  in 
5.1.  A  screening  experiment  did  not  reveal  substantial  differences 
between  the  creep  rates  of  specimens  cut  from  different  panels. 

The  sequence  of  elevated  temperature  creep  and  creep-recovery 
tests,  indicated  in  Fig.  5.12,  was  applied  at  a  moderate  stress  level 
of  a  =  1600  psi  (11  KPa).  The  inverse  compliance  at  t  =  1  min,  i.e. 
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Time 

Fig.  5.12  Creep  (1  hr)  and  Creep-Recovery  (10  hr)  Test  Sequence  at 
Different  Temperatures. 

aQ/e  (1  min)  is  plotted  in  Fig.  5.13.  Moduli  obtained  from  Instron 
tests  are  plotted  on  the  same  graph.* 

The  creep  and  creep-recovery  curves  obtained  during  the  test 
sequence  given  in  Fig.  5.12  are  very  accurately  modelled  by  the  power 
law  equation, 

D(T)  =  Dq(T)  +  D1  tn  (2.19) 

and  the  result  is  shown  in  Fig.  5.14. 


*The  Instron  test  results  were  obtained  using  a  2-in  extensometer , 
whereas  the  creep  and  creep-recovery  results  were  obtained  using 
Micro-Measurements  WK-06-1 25AD-350  strain  gages. 


66 


Temperature  (  deg  C  ) 


50  100  150 


100  200  300 

Temperature  ( deg  F  ) 

Fig.  5.13  Neat  Resin  Modulus  Versus  Temperature. 

The  temperature  dependence  of  DQ(T),  (Fig.  5.15),  indicates  a 
thermorheologically-complex  material.  Figure  5.15  also  shows  that  the 
simple  power  law  [Eq.  (5.21)]  breaks  down  above  300°F  (149°C).  The 
temperature  dependence  of  D-j(T)  versus  temperature  is  shown  in  Fig. 
5.16.  The  irregularity  in  this  curve  identifies  the  w-transition.  The 
transition  is  centered  around  200°F  (93.3°C),  which  is  in  agreement 
with  the  results  obtained  by  Keenan  et  al .  [79],  i.e.  212°F  (100°C),  on 
this  type  of  epoxy  as  well  as  those  reported  by  Griffith  [7]  for  the 
co-transition  of  the  T300/934  composite  system,  220°F  (104°C).  The 
value  obtained  by  Yeow  [5],  140°F  (60°C),  could  be  simply  due  to  the 
combined  effects  of  moisture  absorption  and  uncertainties  in  the  iden¬ 
tification  of  a  small  discontinuity  in  the  thermal  expansion  versus 
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temperature  curve. 

The  master  curve,  obtained  by  horizontal  and  vertical  shifting 
of  the  data  in  Fig.  5.14  is  shown  in  Fig.  5.17  while  the  amount  of 
vertical  versus  horizontal  shift  is  shown  in  Fig.  5.18.  This  latter 
figure  can  be  approximated  by  two  straight  lines  with  slopes  a  and  3 
respectively.  The  discontinuity  in  the  slope  is  an  indication  of  the 
change  from  one  activation  mechanism  to  another.  It  is  well  known  that 
the  horizontal  shift  in  the  glassy  region  can  be  represented  by  an 
Arhenius  type  of  temperature  dependence 


log  ay 


af  n  i  i 

2.303  R  |_T  "  Tr 


where  aF  is  the  activation  energy  (per  mole),  R  is  the  universal  gas 
constant,  and  T^  is  an  arbitrary  reference  temperature.  Figure  5.19 
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Fig.  5.19  Horizontal  Shift  Versus  the  Reciprocal  of  Temperature. 

reveals  a  discontinuity  in  the  log  a^  versus  T"^  curve.  This  again  is 
an  indication  of  a  changing  mechanism.  The  activation  energy  below 
220°C  was  found  to  be  4.76  kcal/gm-mole  while  above  220°F  (104°C)  a 
value  of  23.50  kcal/gm-mole  was  found. 

Instron  stress-strain  curves  (Fig.  5.20)  revealed  a  brittle- 
ductile  transition  in  the  temperature  range  200-300°F*  (93-149°C). 

All  curves  have  an  initially  linear  stress-strain  behavior.  The 
values  of  stress  and  strain  at  the  "linear-nonlinear  transition  point" 
are  plotted  versus  temperature  in  Fig.  5.21.  Based  on  these  results, 
a  test  program  was  developed.  At  a  given  temperature,  the  same 
specimen  was  forced  through  a  series  of  creep  and  creep  recovery  tests 

*Probably  after  the  w-transition  has  decayed. 
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Fig  5.20  Stress-Strain  Curves  for' 934  Resin  Postcured  at  350  F 
(177  C) 
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Fig.  5.21  Values  of  Stress  (o)  and  Strain  (a)  at  the  Linear- 
Nonlinear  Transition  Point. 


as  indicated  in  Fig.  5.22.  A  fresh  specimen  was  used  for  each  new 
temperature.  Test  results  obtained  at  five  different  temperatures, 
with  at  least  eight  subsequent  creep  and  creep-recovery  tests  at  each 
temperature,  up  till  320°F  (160°C)  did  not  show  evidence  of  nonlinear 
viscoelastic  behavior.  As  an  example,  Fig.  5.23  shows  that  the 
compliance  versus  time  evolution  at  the  highest  temperature,  320°F 
(160°C),  is  stress  independent.  Similar  results  were  obtained  for  the 
lower  temperature  tests. 

When  we  proceeded  through  the  loading-unloading  cycles  (Fig.  5.22) 
we  observed  an  accumulation  of  nonrecoverabl e  strain  (Fig.  5.24).  The 
importance  of  this  permanent  strain  accumulation  process  with  respect 
to  the  viscoelastic  deformation  process  can  be  judged  when  we  replot 


1  hr 

1  hr 

1  hr 

o 


Time 

Fig.  5.22  Test  Sequence  of  Creep  and  Creep  Recovery  Tests  at  a  Single 
Temperature  and  Different  Stress  Levels. 
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Fig.  5.24  with  a  vertical  scale  R,  being  defined  as: 


total  permanent  set  accumulated  after  k  creep 

n  =  _ and  creep-recovery  cycles _  , QQ0/ 

the  time  dependent  deformation  at  stress  0 

level 

This  figure  (5.25)  shows  that  the  permanent  set  accumulation  tends  to 
30%  of  the  transient  deformation  at  the  highest  stress  levels.  The 
nonlinearizing  parameters  g ^  and  aa  became  stress  independent  when  we 
subtracted  the  accumulated  strain  (Figs.  5.26  and  5.27).  This  leads 
to  the  conclusion  that  the  matrix  material  behaves  linear  viscoelastic 
at  lower  stress  levels  and  linear  viscoelastic-plastic  at  higher  stress 
levels. 

(MPa) 


10  20  30  40 


Fig.  5.25  Relative  Importance  of  Permanent  Deformation  Compared  to  the 
Viscoelastic  Deformation  for  Different  Stress  Levels. 
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Fig.  5.27  (a)  Shift  Factor  for  934  Resin  with  Subtraction  of 

Accumulated  Permanent  Strain. 


(MPa) 
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Fig.  5.27  (b)  Shift  Factor  for  934  Resin  Without  Subtraction  of 

Accumulated  Permanent  Strain. 
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Chapter  6 

EXPERIMENTAL  RESULTS  AND  DISCUSSION  FOR  T300/934  COMPOSITE 
6.1  Introduction 

Only  creep  measurements  were  made  in  our  previous  studies  [4-13]. 
The  Schapery  nonlinear  viscoelastic  method  described  earlier  could  not 
be  used  as  it  requires  both  creep  and  creep  recovery  data.  The  present 
investigation  was  undertaken  to  obtain  the  proper  information  on 
T300/934  laminates  such  that  the  Schapery  analysis  could  be  developed. 
We  also  wished  to  compare  the  results  of  our  measurements  with  new 
batches  of  T300/934  from  a  new  source,  Fiberite  (see  Appendix  A  for 
curing  details),  and  earlier  batches  supplied  by  Lockheed-Sunnyval e . 

A  further  objective  was  to  compare  the  nonlinear  viscoelastic  proper¬ 
ties  of  the  T300/934  composite  with  those  of  the  934  resin  given  in 
Chapter  5. 

The  plane  stress  constitutive  properties  of  a  linear  orthotropic 
material  is  completely  characterized  by  four  independent  properties. 
These  may  be  the  compliances,  D. ^ ,  D,„  =  D^-j ,  D^  and  Dgg.  The  first 
two,  D-j  i  and  D-^,  can  be  determined  from  a  tensile  test  of  a  uni¬ 
directional  laminate  tensile  coupon  with  the  fibers  in  the  load  direc¬ 
tion  (see  Fig.  6.1).  D^  and  Dgg  can  be  determined  from  tensile  tests 
of  a  unidirectional  laminate  with  the  fibers  transverse  to  the  load 
axis  and  the  fibers  at  10°  to  the  load  axis  respectively.  Performing 
creep  tests  to  express  these  parameters  as  a  function  of  time  permits 
the  calculation  of  orthotropic  linear  viscoelastic  response  for  an 
arbitrary  direction  of  the  fibers. 
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Fig.  6.1  Off-Axis  Test  Geometry  and  Internal  Stress  State. 

Normally,  because  the  graphite  fibers  of  a  G/E  laminate  are  very 
stiff  and  nearly  time  independent,  creep  tests  to  determine  time 
dependent  responses  are  only  performed  for  the  90°  and  10°  cases 
respecti vely. 

If  an  orthotropic  material  is  nonlinear,  D22  is  not  only  a  func¬ 
tion  of  a 2  but  may  be  a  function  of  a-j  and  as  well  [15].  This  gives 

rise  to  an  interaction  among  stress  components  not  present  for  a 
linearly  viscoelastic  orthotropic  material.  In  order  to  account  for 
this  interaction  effect,  various  methods  have  been  employed  usually 
based  on  the  introduction  of  a  flow  rule.  Such  an  approach  is  not 
always  straightforward  due  to  the  possibility  of  various  stress  inter¬ 
actions  at  the  microlevel.  A  simple  method  to  account  for  interaction 
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effects  is  to  replace  the  one-dimensional  stress  by  the  octahedral 
shear  stress.  This  stems  from  the  intuitive  hypothesis  that  creep  is 
caused  by  the  shearing  of  the  polymer  molecules  past  one  another. 

A  somewhat  more  intricate  approach  was  proposed  for  fiber  rein¬ 
forced  materials  by  Lou  and  Schapery  [49]  and  was  in  turn  used  by 
Dillard  [9].  Their  procedure  utilized  the  octahedral  shear  stress  only 
in  the  matrix  phase  as  a  nonlinearizing  parameter.  The  relationship 
between  the  octahedral  stress  is  the  matrix  and  the  plane  stress 

components  (a-pc^T-^)  shown  in  Fig.  6.1  was  given  through  a  micro¬ 
mechanics  model.  The  compliance  model  proposed  by  Dillard  [9]  was 
given  as. 
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An  alternate  approach  would  be  the  use  of  the  free  energy  expan¬ 
sion  of  Eq.  (3.18)  which  could  accommodate  higher  orders  of  the  strain 
or  stress  components  or  their  invariants.  Hahn  and  Tsai  [83]  used  a 
fourth  order  strain  expansion  of  the  complementary  energy  density  to 
handle  shear  nonlinearities.  In  other  words,  temperature  was  excluded 
from  the  Gibbs  free  energy  of  the  system.  We  feel  this  approach  should 
be  extended  and  explored  further. 

Nonlinear  viscoelastic  effect  was  not  observed  in  the  transverse 
direction  for  our  T300/934  material.  During  the  final  compilation  of 
the  results  it  was  realized  that  this  result  should  be  reexamined, 
since  there  is  some  indication  that  the  nonlinear  mechanism  could 
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operate  on  such  a  long  time  scaTe  that  it  becomes  impossible  to 
discriminate  between  materials  response  and  experimental  error  on  the 
basis  of  short  time  accelerated  characterization  results. 

6.2  Results  and  Discussion  on  [90]g$ 

The  laminate  used  in  this  investigation  consisted  of  sixteen  uni¬ 
directional  plies,  with  about  sixty  percent  volume  of  T300  graphite 
fibers  in  a  matrix  of  Fiberite-934  epoxy  resin. 

Tensile  specimens  were  cut  from  a  unidirectional  plate  by  means 
of  a  diamond  cutting  wheel,  postcured  at  350°F  (177°C)  for  4  hours  and 
dessicated  as  discussed  earlier.  A  series  of  elevated  temperature 
creep  and  creep-recovery  tests  were  performed  in  the  same  manner  as 
those  for  the  resin  shown  in  Fig.  5.12.  The  strain  response  was 
measured  using  back  to  back  strai n  gauges  to  avoid  the  effects  of  bend¬ 
ing.  Parenthetically,  it  is  easily  shown  that  an  eccentricity  (e)  of 
1%  of  the  depth  (d)  of  the  rectangular  cross  section  as  shown  in  Fig. 
6.2  causes  a  bending  stress  at  the  gauge  location  which  is  6%  of  the 
average  stress  c^.  (This  is  discussed  in  more  detail  in  Appendix  C.) 


Fig.  6.2  Possible  Bending  of  Tensile  Specimen 
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The  eccentricity  in  creep  tension  tests  has  been  estimated  by 
Isaksson  [84]  to  be  of  the  order  of  5 %  of  the  depth  (d)  causing  a  30% 
deviation  of  the  average  stress  o^. 

The  change  of  the  gauge  factor  with  temperature  was  neglected 
because  it  was  typical  only  1.5%  for  a  temperature  change  from  room 
temperature  up  to  320°F  (177°C). 

Another  experimental  necessity  was  to  avoid  heat  leakage  through 
the  lower  oven  load  train  entrance  as  illustrated  in  Fig.  6.3.  Efforts 
to  insulate  this  entrance  in  the  same  way  as  the  top  load  train 
entrance  caused  interference  with  the  recovery  data.  The  rod  that  was 
used  to  connect  and  disconnect  the  specimen  with  the  lower  and  movable 
end  of  the  load  train  needed  to  be  completely  free.  A  satisfactory 
solution  for  this  problem  was  to  place  a  heavy  mass  of  metal  at  the 
lower  entrance  as  shown  in  Fig.  6.3b.  This  "heat  sink"  contained  an 


Fig.  6.3  Loading  and  Oven  Insulation  Details 
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opening  trough  which  the  connecting  rod  could  move  freely,  but  at  the 
same  time  allowed  a  uniform  oven  temperature  to  be  maintained. 

The  transverse  compliance  as  a  function  of  time  for  different 
temperatures  in  the  range  200-320°F  (93-1 60°C)  is  shown  in  Fig.  6.4. 

An  expanded  vertical  scale  was  used  as  compared  with  the  similar  resin- 
compliance  plot  of  Fig.  5.3,  because  the  continuous  interruptions  of 
the  matrix  by  time  independent  filaments  prohibited  a  large  transverse 
creep  strain  buildup. 

Additional  tests  at  different  stress  levels  were  performed  for  the 
transverse  compliance  as  a  function  of  temperature  as  shown  in  Fig.  6.4. 
However,  these  are  not  shown  as  the  compliance  appeared  to  be  inde¬ 
pendent  of  stress  level,  i.e.,  was  linear,  for  the  time  scale  of  our 
tests  as  will  be  discussed  subsequently. 


Log  Time  (min) 

Fig.  6.4  Transverse  or  D22  Compliance  of  T300/934  as  a  Function  of 
Time  and  Temperature 


83 


The  inverse  compliance  at  t  h  1  min  is  plotted  versus  temperature 
in  Fig.  6.5.  Similar  results  obtained  by  Yeow  [5]  and  Griffith  [7]  are 
plotted  on  the  same  graph,  together  with  a  micromechanics  prediction, 
based  on  the  previously  measured  resin  modulus  as  a  function  of 
temperature. 

Temperature  (degC) 


50  100  150 


100  200  300 

Temperature  (  deg  F  ) 

Fig.  6.5  T300/934  1 /D22  (1  min)  Versus  Temperature 


The  micromechanics  equation  used  to  produce  Fig.  6.5  was  the 
Halpin-Tsai  equation  [54] 

1  +  C*n*v^ 

E2(t,T)  =  Em(t,T)  1  -  n-vf  (6'2 

where  E./E  -  1 

f  m 

°  "  Ef/Em  +  5 

in  which  5  is  a  measure  of  rei nforcement  and  depends  on  the  fiber 
geometry,  packing  geometry  and  loading  conditions  as  discussed  in  [54], 


84 


Equation  6.2  is  valid  only  for  an  elastic  material  but  can  be  extended 
to  the  viscoelastic  problem  using  the  so-called  "quasi-elastic" 
approach.  It  means  that  if  the  solutions  to  the  elastic  field  equations 
are  found,  an  approximate  viscoelastic  solution  can  be  obtained  by  re¬ 
placing  elastic  constants  by  time  dependent  moduli.  Sims  [85]  showed 
that  for  a  creep  time  of  several  hours,  this  quasi -elastic  approach 
deviates  little  from  the  more  theoretically  justified  Laplace  transform 
inversion  approach.  In  Fig.  6.5  it  was  also  tacitly  assumed  that 
E(t)  =  1/D(t),  which,  according  to  Aklonis  and  Tobolsky  [86],  does  not 
introduce  an  error  larger  than  our  experimental  uncertainty  in  deter¬ 
mining  the  compliance. 

As  may  be  observed,  good  agreement  was  obtained  between  the  Hal  pi n- 
Tsai  micromechanics  equation  and  current  and  previous  experimental  work. 
Such  is  the  case  even  though  the  model  is  based  on  isotropic  properties 
on  both  fibers  and  matrix  even  though  the  graphite  fibers  are  known 
to  be  very  anisotropic.  This  confirms  a  statement  by  Whitney  [87] 
which  indicates  that  a  straightforward  substitution  of  appropriate 
isotropic  constants  appears  to  be  adequate  to  account  for  transversely 
isotropic  fiber  properties. 

Micromechanics  solutions  based  upon  an  exact  or  approximate  solu¬ 
tion  to  the  field  equations  have  been  developed,  however,  which  include 
complete  anisotropy  of  the  fibers.  These  include  the  upper  and  lower 
bound  method  of  Bulavs  et  al .  [90]  and  Hashin  [91].  Experimental 
studies  have  been  conducted  by  Ishikawa  et  al .  [88]  and  Thompson  [89] 
which  show  the  Thornel  300  high  strength  fibers  to  be  greatly  aniso¬ 
tropic.  These  results  are  shown  in  Table  6.1. 
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Table  6.1 


I s hi kawa  et  al .  [88] 
(T-300A) 

Thompson  [89] 
(T-300,  wyp  15) 

E11 

32.5 

(MSI) 

32  ±  3  (MSI) 

E22 

2.93 

(MSI) 

/ 

V1 2 

.2 

/ 

v23 

.42 

/ 

G1 2 

6.82 

(MSI) 

/ 

A  comparison  of  Halpin-Tsai  equations  (6.2),  Hashin's  upper  and 
lower  bound  method  with  anisotropic  properties,  Hashin's  equations  with 
equivalent  longitudinal  (E^-, )  and  transverse  (E^t)  fiber  properties, 
Bulav's  square  array  method  and  our  current  experiments  are  given  in 
Fig.  6.6.  As  may  be  observed,  Whitney's  conclusion  is  once  again  con- 
fi rmed. 

Hashin's  upper  and  lower  bound  results  are  shown  in  Fig.  6.7  as 
a  function  of  the  ratio  of  longitudinal  to  transverse  fiber  moduli. 

Equation  6.1  was  modified  to  include  the  micromechanics  approach 
given  by  Eq .  6.2.  Using  the  quasi -elastic  approach  mentioned  pre¬ 
viously  predictions  of  the  time  dependent  transverse  modulus  were  made 
and  these  are  given  together  with  experimental  results  in  Fig.  6.8. 

Good  agreement  was  obtained  which  would  indicate  a  simple  time  dependent 
micromechanics  analysis  to  be  a  reasonable  approach  to  represent  our 
data  over  the  time  range  involved. 
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Fig.  6.8  Comparison  Between  Experimental  (-)  and  Micromechanics  (a) 
Prediction  of  E22  Versus  Time  (T300/934 


An  attempt  was  made  to  model  the  creep  response  with  a  Schapery 
type  power  law  equation 


a(T)_(T) 

g  T  >  0  +  — - — —  n 

ao  0  n  1 

d. 

a 


(6.3) 


where  all  quantities  are  as  defined  previously  in  Equation  3.49.  The 
upper  index  T  indicates  that  the  nonlinearizing  functions  are  evaluated 


in  the  transverse  fiber  direction. 


A  series  of  40  (one  hour)  creep  and  (ten  hour)  creep  recovery  tests 
at  eight  different  stress  levels  and  five  different  temperatures,  using 
a  fresh  specimen  at  each  new  temperature,  did  not  show  any  dependence  of 
the  transverse  compliance  on  the  stress  level.  In  other  words. 


which  means  that  the  transverse  response  was  linear  within  our  time 
frame  and  within  experimental  error. 
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The  stress-strain  diagram  (Fig.  6.9)  obtained  from  an  Instron 
test  at  100°F  (37°C)  also  shows  little  dissipation  and  a  fracture  strain 
of  about  5000  microstrain. 


Strain  (%) 

Fig.  6.9  Stress-Strain  Curve  (Solid  Line)  for  [90]gs  T300/934  Laminate 
at  Room  Temperature  (0.05  in/min)  Compared  with  Linear 
Response  (Dashed  Line). 

Equation  (6.3)  reduces  to 

D(t)  =  Dq  +  D1  tn  (6.4) 

Figure  6.10  shows  the  temperature  dependence  of  the  first  term  (D  )  in 
Eq.  (6.4).  It  is  clear  that  the  dependence  is  almost  perfectly  linear 
in  the  temperature  interval  which  we  studied.  The  temperature  dependence 
of  the  coefficient  D-|  has  a  pronounced  nonlinear  variation  as  can  be 
seen  in  Fig.  6.11.  Notice  that  the  w-mechanism  which  was  identified  in 
the  neat  resin  (Fig.  5.2)  is  also  clearly  visible  in  Fig.  6.11.  However 
the  minimum  in  the  curve  is  now  approximately  20°F  (=  10°C)  higher. 
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Such  differences  have  been  reported  for  the  primary  transition  by  a 
number  of  investigators  [76,79,92],  but  the  cause  is  unclear. 

The  long  time  creep  data  plotted  in  Fig.  6.12  reveals  that  our 
experimental  time  scale  for  short  time  creep  (1  hour)  was  probably  too 
short  to  be  able  to  discriminate  between  stress  dependence  of  the 
compliance  curve  and  small  experimental  irregularities  such  as  mis¬ 
alignment  in  the  load  train,  fluctuation  in  temperature,  etc.  Or,  in 
other  words,  the  shorter  the  test  period,  the  more  precise  the  test 
conditions  need  be  in  order  to  allow  stress  dependent  discrimination  of 
the  compliance  curve.  Thus  our  earlier  conclusion  about  the  linear 
viscoelastic  or  stress  independence  of  our  transverse  compliance  data 
for  T300/934  needs  to  be  reexamined  over  a  longer  time  scale.  However 


Fig.  6.12  Comparison  of  Transverse  Compliance  Versus  Time  for  Moderate 
(2.76  Ksi)  19  Mpa)  and  High  (5.18  Ksi  =  35.7  Mpa)  Stresses 
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for  a  first  approximation,  linear  would  be  a  reasonable  assumption. 

6.3  Results  and  Discussion  on  [10]gs 
6.3.1  Experimental  Background 

An  efficient  intralaminar  (in  plane)  shear  characterization  tech¬ 
nique  known  as  the  10-deg  off-axis  test  has  been  proposed  and  analyzed 
by  C.  C.  Chamis  and  J.  H.  Sinclair  [93].  Their  analysis  shows  that  the 
major  contribution  to  fracture  comes  from  intralaminar  shear,  i.e.  the 
configuration  activates  a  long  continuous  shear  path,  along  which  load 
diffusion  into  the  reinforcing  fibers  occurs.  One  difficulty  with  the 
off-axis  test  is  that  shear  coupling  effects  are  a  function  of  the  load 
to  fiber  angle  or  the  aspect  ratio.  The  Pagano-Hal pin  [94]  shear 
coupling  analysis  was  used  to  generate  Fig.  6.13.  Their  analysis  sug¬ 
gest  that  for  an  aspect  ratio  or  off-axis  angle  of  10°,  the  shear 
modulus,  Gi2^approx  ’  is  aPProximately  10°^  higher  than  the  true  shear 
modulus  G-|2-  This  is  confirmed  by  the  inverse  shear  compliance  versus 
temperature  plot  of  Fig.  6.14  which  indicates  a  measured  value  some¬ 
what  above  .8  MSI  at  room  temperature,  while  a  widely  published  value 
is  closer  to  .7  MSI. 

Two  back-to-back  45  deg  rosette  strain  gages  [WK-06-062R-350] 
were  used  to  eliminate  out  of  plane  bending  from  test  results.  The 
orientation  of  the  gages  with  respect  to  the  fiber  axis  was  checked 
under  the  microscope.  Misalignments  were  taken  into  account  during 
subsequent  data  reduction.  Figure  6.15  shows  the  orientation  of  the 
strain  gage  rosette  as  well  as  the  manner  in  which  aluminum  tabs  were 
adhesively  bonded  to  each  specimen. 
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Shear  Modulus  Ratio,  G|2  / G|2  )appar 


Fig.  6.15  Specimen,  Strain  Gage  and  Tab  Configuration. 


Each  specimen  was  postcured  with  the  same  temperature  cycle  as 
were  the  transverse  and  neat  resin  coupons. 

The  shear  compliance  Dg^(t)  for  a  10°  specimen  loaded  in  tension 
as  illustrated  schematically  in  Fig.  6.1  is  given  by  the  following 


expression 


D66(t) 


w-h.(-.598  ex(t)  +  1.879  e^t) 

.171  a.. 


-  1.282  eft)) 


where  ex(t),  e^(t)  and  Sy(t)  are  the  strains  measured  along  the 
longitudinal,  45°  and  transverse  directions  of  the  specimen,  respectively. 
The  angle  between  the  measured  strain  directions  and  the  fiber-direction 


of  the  material  are  10°,  45°  and  100°  respectively.  W  and  h  represent 
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the  width  and  thickness  of  the  tensile  coupon.  The  strains  measured 
along  the  x  and  45°  direction  were  positive,  while  the  strain  measured 
along  the  y-direction  was  negative.  It  should  be  noted,  of  course,  that 
if  the  rectangular  rosette  is  mounted  with  the  90°  and  45°  arms  dif¬ 
ferent  than  shown  in  Fig.  6.15,  different  signs  will  be  obtained  and 
the  preceding  equation  for  Dg6  must  be  altered  appropriately.  The  shear 
compliance  Dg6  is  thus  easily  obtained  through  simultaneous  measurement 

of  e  (t),  6/,c(t)  and  e  (t)  but  caution  should  be  taken  with  respect  to 
x  o  y 

the  gage  orientation. 

6.3.2  Experimental  Results 

The  shear  compliance  versus  log  t  curves  obtained  at  four  dif¬ 
ferent  temperatures  are  shown  in  Figs.  6.16a,  b,  c  and  d.  They  clearly 
indicate  the  behavior  to  be  dependent  upon  stress.  In  the  context  of 
nonlinear  viscoelasticity,  as  discussed  earlier  in  Chapter  3,  this  means 
that  a  set  of  molecular  processes  which  are  sensitive  to  stress  must 
also  contribute  to  the  observed  deformation  in  addition  to  those  nor¬ 
mally  associated  with  linear  viscoelasticity. 

A  definite  accumulation  of  nonrecoverable  permanent  strain  was 
observed  during  the  creep  and  creep-recovery  test  sequences.  The  amount 
obtained  at  246°F  ( 1 1 9° C )  is  shown  as  an  example  in  Fig.  6.17.  Thus  the 
molecular  structure  moves  out  of  the  positions  which  it  occupied  on 
release  of  the  stress  but  does  not  reach  its  former  spatial  arrangement 
under  the  action  of  the  restoring  forces.  One  possibility  is  that  this 
former  spatial  arrangement  is  unattainable  due  to  a  "plasticity"  type 
structural  change  in  the  matrix  phase.  Thus,  for  our  material  it  seems 
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Log  Time  ( min  ) 


Fig.  6.16a  Shear  Compliance  Versus  Time  at  246°F  (119°C). 


Log  Time  (  min  ) 


Fig.  6.16b  Shear  Compliance  Versus  Time  at  300°F  (149°  C). 
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Fig.  6.16c  Shear  Compliance  at  320°F  (160°C). 


o  0.0  1.0  2.0  3.0  4.0 


Log  Time  (min) 

Fig.  6.16d  Shear  Compliance  Versus  Time  at  335°F  (168°C) 
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6.0  12.0  18.0  24.0  30.0  36.0 


Fig.  6.17  Permanent  Strain  Accumulation  During  10  Deg.  Off  Axis  Test 
at  246°F  (11 9°C) . 

not  possible  for  nonlinear  viscoelasticity  to  occur  without  involving 
other  irreversible  deformation  mechanisms.  This  conclusion  disagrees 
with  results  obtained  by  Ferry  [28]  on  PMMA.  Shear  creep  deformation 
up  to  5%  were  fully  recoverable  although  the  linear  viscoelastic  limit 
was  below  1%. 

A  number  of  investigators  mechanically  condition  specimen  prior 
to  creep  testing.  The  reason  is  that  during  the  first  creep  test  a 
certain  amount  of  deformation  takes  place  which  does  not  recover  upon 
removal  of  the  stress.  In  subsequent  creep  tests  on  the  same  specimen 
at  the  same  or  smaller  load  there  is  no  such  nonrecoverabl e  creep.  The 
stress-strain  relationship  thus  becomes  repeatable.  The  material  is 
therefore  considered  "mechanically  conditioned."  Leaderman  [22]  showed 
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that  the  permanent  set  in  fibers  disappears  when  the  fiber  was  steamed 
under  no  load  and  then  dried.  Thus  no  real  structural  change  or 
plasticity  was  introduced  by  the  conditioning  process. 

This  author  fears,  however,  that  a  significant  volume  of  the 
matrix  material  might  be  yielded  by  the  very  inhomogeneous  stress  fields 
caused  by  the  geometrical  nature  of  the  composite  material.  Evidence 
of  such  a  matrix  yielding  phenomenon  in  shear  deformation,  based  on  a 
finite  element  analysis,  was  published  by  Foye  [95].  His  results  shown 
in  Fig.  6.18a  give  some  indication  of  the  progress  of  the  plastic  zone 
boundaries  due  to  shear  loading  in  repeating  elements  of  a  graphite- 
epoxy  composite  with  50%  fiber  by  volume.  Yielding  initiates  at  the 
interface  in  an  area  enclosed  by  plastic  zone  boundary  number  1.  The 
plastic  zone  boundaries  numbered  2  through  8  correspond  to  increasing 
postyielding  load  levels.  Figure  6.18b,  on  the  other  hand,  shows  that 
in  transverse  deformation,  yielding  first  occurs  at  a  point  in  the 
matrix  midway  between  two  adjacent  fibers.  Foye  also  showed  that  initial 
yielding  of  the  composite  during  the  transverse  deformation  mode  occurs 
at  slightly  lower  average  normal  stress  levels  than  the  bulk  matrix. 

In  the  shear  deformation  mode  the  composites  yielded  at  much  lower 
average  shear  stress  levels  than  the  bulk  matrix. 

In  conclusion,  preconditioning  was  not  made  part  of  the  testing 
program  since  it  was  feared  that  the  stress-strain  relationship  would 
become  repeatable  at  the  expense  of  changing  the  matrix  phase  to  an 
entirely  new  material  in  the  sense  of  a  plasticity  type  shakedown  model. 


Matrix 


Fiber 


Fig.  6.18a  Inhomogeneous  Stress  Distribution  Due  to  Shear  Deformation 
Due  to  Shear  Deformation  After  Foye  [95]. 
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Fiber 


Fig.  6.18b  Inhomogeneous  Stress  Distribution  Due  to  Transverse  Deforma¬ 
tion  After  Foye  [95]. 


The  creep  response  in  the  shear  deformation  mode  was  modeled  with 


ioo 


the  power  law 


6.3.3  Time-Stress  Superposition  Principle  (TSSP) 

A  very  useful  viscoelastic  concept  is  that  time  and  stress  are 
equivalent  for  describing  viscoelastic  behavior  [7].  This  is  due  to  the 
fact  that  an  increased  tensile  stress  shifts  the  creep  spectra  towards 
shorter  times  as  shown  in  Fig.  6.20a.  Thus,  as  shown  in  Fig.  6.20b, 
creep  observed  for  short  times  at  a  given  stress  x-j  is  identical  with 
creep  observed  for  longer  times  at  a  lower  stress  xQ  provided  that  these 
data  are  shifted  on  a  logarithmic  time  axis.  Similarly,  portions  of  the 

creep  curve  for  x  =  xQ  can  be  observed  for  stresses  x  =  Tptg .  as 

indicated  by  the  curved  segments  shown  in  Fig.  6.20b.  These  curved 
segments  can  then  be  shifted  along  the  log  time  axis  to  construct  a 
composite  curve  or  the  sigmoidal  master  curve  applicable  for  the  given 
stress  level  r  ,  extending  over  many  decades  of  time. 
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Fig.  6.19c  Nonlinear  Parameter  92(1). 
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Fig.  6.1 9d  Nonlinear  Parameter  a^t). 
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Fig.  6.20b  Dependence  of  Compliance  D  on  log  t  and  Stress  t. 


The  logarithm  of  the  shift  factor  aT  for  a  particular  curved 
segment  is  the  horizontal  shift  necessary  to  allow  it  to  join  smoothly 
into  the  master  curve.  Each  of  the  curved  segments  in  Fig.  6.20b  can 
be  described  by 
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This  same  equation  has  a  limited  capability  to  represent  a  master 
curve  using  the  form. 
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(6.8) 


The  reason  for  this  limitation  is  due  to  the  approximate  nature  of  the 
power  law  used  and  was  discussed  earlier  in  Chapter  2. 

The  time-stress-superposition  principle  (TSSP)  described  above 
was  used  to  convert  the  data  of  Figs.  6.16a,  b,  c  and  d  to  master 
curves  for  each  temperature  and  are  given  in  Figs.  6.21a,  b,  c  and  d. 
Also,  equation  (6.8)  was  used  to  analytically  determine  each  master 
curve  and  the  resulting  shift  parameters.  The  symbols  represent  the 
master  curves  produced  graphically  while  the  solid  line  represents  the 
analytically  produced  master  curve.  Comparison  of  the  shift  factors 
found  by  the  two  methods  are  shown  in  Figs.  6.22a,  b  and  c.  A  three- 
dimensional  version  of  these  results  is  shown  in  Fig.  6.23.  As  may  be 
observed,  excellent  agreement  was  obtained  between  graphically  and 
analytically  produced  master  curves. 


105 


c n 
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Fig.  6.21a  Sg6  Master  Curve  for  T300/934  at  246°F  (119°C) 


Fig.  6.21b  Master  Curve  for  T300/934  at  300°F  (149°C) 
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6.22a  Comparison  of  Graphical  and  Analytical  Obtained  Horizontal 
Shift  Factor  .10  deg  -  off  -  axis  for  246°F  (119°C). 
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6.22b  Comparison  of  Graphical  and  Analytical  Obtained  Horizontal 
Shift  Factor  .10  deg  -  off  -  axis  for  300°F  (149°C). 


m 


Shift  Factor 


6.22c  Comparison  of  Graphical  and  Analytical  Obtained  Horizontal 
Shift  Factor  .  10  deg  -  off  -  axis  for  335°F  (168°C). 
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6.23  Graphical  and  Computed  Shift  Function  Surface  for  T300/934. 


Stress  Activated  Shift 


Often,  a  vertical  shift  as  well  as  a  horizontal  shift  is  neces¬ 
sary  to  produce  a  smooth  master  curve.  In  fact,  even  a  rotation  of 
individual  compliance  curves  may  be  necessary.  The  type  of  shifting 
needed  depends  on  the  behavior  of  the  limiting  compliances  Du  and  DR 
as  a  function  of  stress  as  shown  schematically  in  Fig.  6.24.  A 

vertical  shift  of  DT  -  D  T°  produces  the  dotted  durve  Dx(t)  which  can 

u  u  r  p 

be  superimposed  on  DT°(t). 

A  rotation  would  be  necessary  when  the  limiting  compliances  are 
affected  differently  by  the  stress  intensity,  i.e.,  when 

V  -  D  T 

—7 - -  i  1  (6.9) 

V°  -  °u  ° 


Stress  Log  Time  (min) 


Fig.  6.24  Illustration  of  the  Stress-Shi fting-Procedure  (Based  on 
McCrum  et  al .  [96] ) . 


In  our  case,  a  small  vertical  shift  was  found  necessary  at  246°F.  This 
shift  however  was  nonexistent  at  higher  temperatures  as  shown  in  Figs. 
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6.25a,  b  and  c,  where  the  vertical  shift  is  plotted  versus  the  hori¬ 
zontal  shift. 

It  was  thought  initially  that  the  true  vertical  shift 
(Dut  -  DuT°)>  given  in  Fig.  6.24,  was  equivalent  to  (gQ  -  1)DQ.  This 
however  leads  to  an  overprediction  due  to  the  fact  that  the  sigmoidal 
compliance  curves  change  location  with  increasing  stress.  The  impli¬ 
cation  is  that  even  when  the  limiting  compliance  values  are  completely 
insensitive  to  stress  (which  turns  out  to  be  the  case  in  this  investi¬ 
gation)  we  still  would  record  a  gQ  value  different  from  unity  since  the 
retardation  spectrum  involves  rapid  creep  on  the  time  scale  of  micro¬ 
seconds.  In  conclusion,  very  accurate  short  time  limiting  compliance 
data  is  needed  to  obtain  a  proper  estimate  of  the  true  vertical  shift. 

As  a  final  conclusion,  a  long  time  control  test  remains  the  only 
criterion  for  checking  a  particular  accelerated  test  method.  These 
long  term  creep  tests  are  currently  being  carried  out  and  will  dictate 
the  adaptations  necessary  to  achieve  a  proper  mathematical  modeling  for 
the  prediction  of  long  term  properties  from  short  term  test  results. 
These  final  modeling  equations  should  in  some  way  reflect  the  mechanism 
of  inelastic  phenomena  operating  at  the  microstructural  level. 
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Chapter  7 

EXPERIMENTAL  RESULTS  AND  DISCUSSION  ON  DELAYED  FAILURE 
AND  LAMINATE  RESPONSE 


7.1  Introduction 

The  usual  concept  of  elastic  deviatoric  strain  energy  in  which 
the  Tsai -Hi  11  anisotropic  failure  criterion  is  based  cannot  be  carried 
over  directly  to  viscoelastic  materials  without  modifications  to 
account  for  dissipation  mechanisms.  However,  a  better  method  to  in¬ 
clude  dissipation  as  well  as  the  entire  history  of  deformation  is  the 
concept  of  free  energy  which  can  be  used  to  define  a  "critical"  stress 
state  at  the  point  of  failure.  A  delayed  failure  model  based  upon 
these  concepts  was  developed  in  Chapter  4. 

The  delayed  failure  model  requires  information  on  both  creep 
and  on  creep  rupture.  An  existing  creep  rupture  data  base  for  [90 ]gs 
and  [60]gs  orientations  of  T300/934  was  available  from  the  work  of 
Griffith  [7].  [90]gs  creep  data  is  given  in  Chapter  6.  The  needed 

[60]gs  creep  response  to  compliment  the  [60]gs  rupture  data  can  be 
inferred  from  the  [90]gs  and  [10]gs  data  given  in  Chapter  6.  However, 
it  was  decided  to  obtain  the  needed  [60]gs  creep  response  from  new 
specimen  because  inferred  response  would  introduce  a  hypothetical 
octahedral  stress  interaction  factor,  t. 

7.2  Delayed  Failure 

7.2.1  Analysis  on  [90]gs  Results 

The  compliance  versus  time  curve  at  320°F  (160°C)  and  1580  psi 
(10.88  MPa)  (Fig.  6.8)  was  fitted  with  a  six-element  Kelvin  model, 
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using  a  collocation  scheme  proposed  by  Schapery  [97].  The  resulting 
compliance  D(t)  can  be  expressed  as: 

D(t)  =  .71  +  .041  [1  -  exp  (-t/.01)]  + 

+  .041  [1  -  exp  ( -t/1 ) ]  +  .060  [1  -  exp  ( - t/1 0 ) ] 

+  .049  [1  -  exp  ( -t/1 00 ) ]  +  .164  [1  -  exp  ( -t/1 000 ) ] 

+  .158  [1  -  exp  (-t/10000)]  x  10"6  psi"1  (7.1) 

Based  on  this  identification  scheme,  we  can  use  Eqs.  (4.9  -  4.11) 
to  calculate  the  total  creep  energy,  stored  free  energy  and  dissipated 
energy  as  a  function  of  time  as  shown  in  Fig.  7.1.  A  gradual  increase 
of  the  total  dissipation  with  time  originates  from  the  energy  con¬ 
sumption  in  the  subsequent  viscous  elements.  This  dissipation  has  to 
be  subtracted  from  the  total  creep  energy,  which  is  in  this  case  equal 
to  the  potential  energy  loss  of  the  external  loading,  in  order  to  ob¬ 
tain  the  total  stored  free  energy. 

An  identical  evolution  of  the  free  energy  was  also  obtained  based 
on  the  power  law  parameters  (D  ,D.j,n)  which  made  it  possible  to 
evaluate  Eq.  (4.17)  as  a  function  of  time.  This  procedure  has  a 
definite  advantage  in  that  it  avoids  a  collocation  scheme. 

A  comparison  of  the  experimental  creep-rupture  results  and  the 
free  energy  based  prediction  shows  that  very  good  agreement  was  ob¬ 
tained  as  shown  in  Fig.  7.2.  The  full  line  (-»— )  represents  a  given 
constant  free  energy  line.  Thus  a  creep  rupture  result  is  obtained 
from  creep  compliance  data  at  the  expense  of  only  one  additional 
parameter  (w).  Earlier  attempts,  in  which  we  tried  to  correlate  the 
total  creep-energy  with  creep  to  rupture  data,  gave  very  poor  agreement. 


114 


Transverse  Stress  ( Ksi )  ^  Energy  /  ( applied  stress )**  2 


1.0  2.0  3.0  4.0 


Log  Time  (min) 


T300/934:  Evolution  of  Total  Creep  Energy,  Total  Stored 
Energy  and  Dissipated  Energy  as  a  Function  of  Time  for 
[ 90 ] 8s  at  320  deg  F  (160°  C). 


Log  Time  (min) 


Fig.  7.2  T300/934:  Comparison  of  Free  Energy  Based  Creep  to  Rupture 
Prediction  with  Experimental  Results  Obtained  on  [90]R 
at  320  deg  F. 


115 


7.2.2  Analysis  of  [60]gs  Results 

Long  time  creep  and  creep  recovery  data  were  obtained  on  a  60-deg 
off  axis  specimen.  The  details  on  postcuring,  strain  gaging  and  test¬ 
ing  were  identical  with  those  already  discussed  in  Chapters  5  and  6. 

A  plot  of  the  compliance  change  with  log  t  at  a  stress  level  of 
72%  of  ultimate  is  shown  in  Fig.  7.3.  A  result  obtained  by  Griffith 
[7]  on  the  same  orientation  but  at  55%  of  ultimate  is  also  shown  in 
this  figure.  The  curvature  of  these  plots  is  clearly  affected  by  the 
stress  intensity.  In  addition,  however,  a  rather  large  vertical  shift 
would  be  needed  to  superimpose  both  curves.  This  could  be  due  to 
matrix  plasticity,  since  a  rather  large  permanent  deformation  remained 
after  creep  recovery  (Fig.  7.4).  Martirosyan  [98]  gave  evidence  of 
such  results  which  are  reproduced  in  Fig.  7.5.  The  nonrecoverable 
creep  deformation  is  seen  to  behave  as  an  increasing-decreasing  func¬ 
tion  with  fiber  angle.  These  bell-shaped  curves  are  parameterized  by 
a  number  that  indicates  the  stress  intensity  as  a  fraction  of  ultimate. 
Obviously  a  large  nonrecovered  strain  was  observed  which  varies  very 
nonli nearly  with  stress  intensity. 

The  creep  response  curve  at  72%  of  ultimate  was  again  modeled 
with  a  six  element  Kelvin  model,  with  an  additional  free  dashpot  in 
series,  as  shown  in  Fig.  7.6. 

The  free  energy  at  rupture  is  seen  to  be  constant  for  delayed 
failures  in  excess  of  30  min.  as  shown  in  Fig.  7.7.  Exactly  the  same 
free  energy  evolution  with  time  was  found  again  through  the  power  law 
approximation.  The  power  law  parameters  were  evaluated,  however, 
after  subtraction  of  the  free  dashpot  contribution. 


116 


Fig.  7.3  Long  Time  Creep  Compliance  of  [60]Q_  T300/934  at  320°F 
(160°C) .  BS 


Fig.  7.4  Typical  Creep  Curve  for  [60]gs  T300/934  at  320°F  (160°C). 
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■7.6  Kelvin  Element  Used  to  Model  Creep  Respon 
for  [60]gs  at  320°F  (160°C). 


Time  ( min ) 

Fig.  7,7  Free  Energy  at  Rupture  for  [60]R_  T300/934 
at  320°F  (1 60°C) . 

A  reasonable  free  energy  based  creep-rupture  prediction  was 
obtained,  as  seen  in  Fig.  7.8.  Final  conclusions  however  cannot  be 
made  yet,  due  to  the  limited  number  of  experimental  data  points  which 
are  currently  available. 

7.3  Creep  Response  on  [90/±45/90]2S ,  Compared  with  Creep  response  of 
[90]8s  and  [±45]4s  at  320°F  (160°cT 

It  can  be  expected,  that  a  [90/±45/90]2s  laminate  configuration 
will  operate  at  a  lower  creep  rate  than  a  [±45]^s»  since  the  latter 
contains  only  two  fiber  orientations  which  allows  scissoring  action 
as  may  be  visualized  in  Fig.  7.9.  The  former  contains  three  fiber 
orientations  and  thus  a  fibre  truss  network  is  formed  as  shown  in 
Fig.  7.9  which  prevents  scissoring. 
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Fig.  7 
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Rupture  Time  (min) 

Fig.  7.8  Free  Energy  at  Rupture  Predictions  for  [60] 
T300/934  at  320°F  (160°C). 


9  Schematic  of  [±45]  and  [90/±45/90]  Laminates  with  Fiber 
Truss  Network  in  Latter. 


Figure  7.10  compares  the  creep  responses  of  three  laminates.  It 
is  qualitatively  clear  that  the  [±45]^s  fiber  orientation  shows  the 
most  transient  creep.  This  can  also  be  quantified  when  we  compare  the 
power  law  coefficients,  D-j ,  which  are  a  function  of  stress,  as  given 
in  Figs.  7.11.  A  comparison  of  the  [90/145/90],^  laminate  results 
with  those  for  the  [90]gs  and  [±45]^s  laminates  reveals  that  the  creep 
rate  for  the  former  is  an  order  of  magnitude  less  than  that  for  the 
other  orientations.  The  initial  compliance  however  decreased  only  by 
35%  as  compared  to  the  [±45]^s  orientation. 


0.0  1 .0  2.0 
Log  Time  (  min ) 

Fig.  7,10  Creep  Compliance  for  Three  T300/934  Laminates. 


In  conclusion,  adding  90' s  to  a  ±  45  orientation  increases  the 
modulus  by  only  35%  but  decreases  creep  sensitivity  by  an  order  of 
magnitude. 

Prior  to  testing  the  [±45 ]4$  laminate,  we  also  put  a  strain  gage 
transverse  to  the  loading  direction.  This  information  allows  us  to 
obtain  the  shear  compliance  [99]  which  is  given  by 


2(6, 
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(7.2) 


It  was  observed,  however,  that  e  (t)  and  e  (t)  had  approximately  the 

x  y 

same  magnitude.  This  led  to  untolerable  large  errors  and  even  in¬ 
consistencies  in  the  creep  recovery  data. 


7.4  Numerical  Prediction  of  [90/±45/90]2S  Laminate  Behavior 

The  incremental  nonlinear  viscoelastic  lamination  theory-program, 
which  was  written  and  documented  by  Dillard  [7  ]  during  an  earlier 
phase  of  our  research  program,  was  used  to  generate  the  creep  response 
of  a  [90/±45/90]2S  laminate. 

A  comparison  of  the  numerical  prediction  with  the  experimental 
result  at  70%  of  ultimate  at  320°F  (160°C)  is  shown,  on  an  expanded 
scale,  in  Fig.  7.12.  Reasonable  agreement  between  predicted  and 
measured  compliance  was  obtained.  Creep  rupture-time,  however,  was 
even  more  seriously  underpredicted,  indicating  a  basic  deficiency  of 
the  numerical  approach  which  lies  in  its  inability  to  reflect  the 
microstructural  effects  of  the  laminate.  This  microstructural  behavior 
was  deduced  using  edge  replicas  [100],  as  shown  in  Fig.  7.13.  An 
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attempt  to  model  the  stress  transfer  mechanism  involved  by  means  of 
a  shear-lag  model  [101],  which  was  simple  enough  for  subsequent 
implementation  in  the  laminate  analysis  program,  turned  out  to  be  un¬ 
successful  . 

Comparison  of  the  creep  response  at  97%  of  ultimate,  also  as 
shown  in  Fig.  7.12,  with  the  response  at  70%  of  ultimate  indicates  a 
large  vertical  shift.  The  edge  replica,  obtained  after  unloading, 
however,  revealed  a  very  different  crack  spacing  between  results  at 
the  two  stress  levels.  The  crack  spacing  in  the  90-deg.  oriented 
lamina  is  considerably  less  than  in  the  previous  case  and  some 
cracking  in  the  45-deg.  oriented  lamina  became  visible  at  the  higher 
load  level . 

Figure  7.14  shows  the  irregularities  in  the  strain  reading  at 
the  beginning  of  the  creep  rupture  experiment  at  97%  of  ultimate. 

This  indicates  a  sequence  of  events  at  the  microstructural  level, 
which  could  be  deterministic  in  the  selection  process  of  the  stress 
transfer  mechanism. 

More  understanding  in  this  area  should  be  generated,  including 
the  answer  on  the  question  of  whether  or  not  crack  healing  sets  in 
under  certain  unloading  or  strain-recovery  conditions. 
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Creep  Time  (mim) 


Fig.  7.14  Creep  Curve  at  97%  of  Ultimate. 
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Chapter  8 

SUMMARY  AND  RECOMMENDATIONS  FOR  FURTHER  STUDY 

Different  constitutive  formulations  able  to  model  material  non- 
linearities  in  viscoelastic  systems  have  been  discussed.  Nonlinear 
viscoelasticity  models  are  often  quite  complicated,  creating  great 
experimental  difficulties  for  the  determination  of  material  parameters. 
It  was  concluded  that  the  approach  advocated  by  Schapery  is  sufficiently 
general  to  account  for  stress-activated  nonlinear  behavior  in  polymers 
and  composites.  Particularly  the  following  experimentally  observed 
phenomena  were  easily  modeled: 

•  An  initial  elastic  creep-recovery  which  is  greater  than  the 
initial  elastic  creep  response 

•  A  recovery  behavior  which  is  overpredicted  by  linear  visco¬ 
elastic  considerations  (i.e.  the  actual  recovery  curve  in  the 
nonlinear  range  is  flatter  than  in  the  linear  range) 

Creep  and  creep-recovery  experiments  which  were  carried  out  on 
the  neat  934-resin  did  not  reveal  any  stress  dependent  change  in 
materials  response  other  than  plasticity.  The  analysis  of  time- 
dependent  response,  which  is  exclusively  connected  with  specifics  of 
the  molecular  structure  revealed  the  existence  of  a  secondary  transi¬ 
tion-temperature  located  between  200  and  250°F  (93-1 21 °C). 

Essentially  the  same  series  of  experiments  were  carried  out  on 
the  T300/934  graphite-epoxy  system.  The  transverse  compliance  D^t) 
was  found  to  be  very  stress  insensitive.  The  shear  compliance  Dgg(t), 
on  the  other  hand,  turned  out  to  be  a  sensitive  function  of  stress. 
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The  role  of  yielding  and  subyielding  of  the  matrix  phase,  due  to  dif¬ 
ferences  in  loading  path  (transverse  loading  vs.  shear  loading)  was 
discussed.  Two  uncoupled  equations  were  fit  to  the  data  D22(t)  and 
Dgg(t)  resPecti vely .  The  calibration  of  the  material  parameters 
(Do,C,n,go,g^ ,g2,a  )  was  done  by  a  computer  routine.  It  is  well  known 
however  that  mathematical  models  fit  data  well  when  used  on  the  stress 
paths  which  were  applied  for  parameter  calibration  purposes.  In  order 
to  account  for  stress  interaction  along  general  loading  paths,  the 
principle  of  octahedral  stress,  was  borrowed  from  nonlinear  approaches 
in  isotropic  viscoelasticity.  The  validity  of  this  approach  however 
needs  to  be  carefully  checked.  The  power  law  creep  exponent  (n),  was 
found  to  be  .23  and  independent  of  stress  and  temperature.  The 
secondary  transition,  which  was  identified  in  the  resin,  was  also 
identified  in  the  composite,  but  at  a  slightly  higher  temperature. 

Short  time  data  could  be  shifted  graphically  into  a  unique  "master 
curve."  The  horizontal  shift  was  found  to  be  in  very  good  agreement 
with  the  analytical  shift  (log  a^).  The  vertical  shift  was  very  small, 
and  became  even  nonexistent  at  higher  temperatures.  It  was  shown  that 
vertical  shifting  is  related  to  the  behavior  of  the  limiting  compliances. 
An  interesting  point  in  this  context  is  that  Urzhumtsev  [102]  presents 
evidence  that  the  time-stress-superposition  principle,  when  used  to 
shift  short  term  compliance  curves  obtained  on  porous  polyurethane  is 
in  excellent  agreement  with  a  long  time  control  experiment  whereas  the 
time-temperature  superposition  leads  to  very  serious  overpredictions 
of  long  time  compliance  data. 
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It  should  be  noted  that  no  mechanical  conditioning  was  applied, 
since  conditioning  should  be  used  to  isolate  delayed  elastic  response, 
which  is  something  that  has  never  been  proven  for  fiber  reinforced 
composites  (at  least  not  to  this  author's  knowledge). 

Creep  compliance  data  were  merged  into  a  free  energy-based 
failure  criterion  and  reasonable  creep-rupture  predictions  were  ob¬ 
tained.  The  creep  response  of  a  [90/±45/90]2s  laminate  was  found  to 
be  an  order  of  magnitude  less  than  the  creep  response  of  [90]gs  and 
[±45]^s,  while  the  modulus  increased  only  by  a  factor  of  1.35.  A 
reasonable  creep  compliance  versus  time  prediction  was  obtained  by  means 
of  a  nonlinear  viscoelastic  laminate  analysis  program.  It  is  felt, 
however,  that  a  damage  model  that  operates  on  the  minimechanics  level 
of  the  problem  would  give  a  better  agreement  than  the  cumulative 
damage  model  which  was  actually  used  [9].  Promising  ideas  in  this 
direction  were  published  recently  by  Nuismer  et  al .  [103].  They 
developed  a  realistic  and  simple  scheme  that  accounts  for  influences 
on  the  constitutive  behavior  of  given  matrix  dominated  laminas  through 
the  lamina-laminate  interaction  after  damage  initiation.  Another,  more 
fundamental,  approach  that  describes  the  influence  of  damage  through  a 
continuum  vector  field  has  been  pursued,  within  a  continuum  thermo¬ 
dynamics  framework  by  Talreja  [104]. 

Based  on  this  current  research  effort,  the  following  recommenda¬ 
tions  for  further  study  seem  appropriate: 

•  The  feasibility  of  mechanical  conditioning,  in  case  of 
composites,  should  be  thoroughly  checked. 
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Test  programs  for  boron-epoxy  by  Cole  and  Pipes  [14]  and  for 
graphite  polymide  by  Pindera  and  Herakovich  [15]  produced  evi¬ 
dence  that  the  octahedral  shear  stress  is  not  a  valid 
representation  of  stress  mode  interaction.  Pindera  and 
Herakovich  showed  that,  in  order  to  get  access  to  the  actual 
stress  interaction,  3  off-axis  tests  were  needed.  An 
identification  of  the  actual  representation  for  stress  mode 
interaction  in  case  of  graphite-epoxy  T300/934  would  at  least 
give  indications  on  its  relative  importance  with  respect  to 
nonlinearities  in  transverse  deformation  and  shear. 

The  free  energy  based  failure  criterion  together  with  the 
Tsai -Hi  11  criterion,  could  be  developed  to  predict  delayed 
off-axis  failure,  based  on  transverse  and  shear  compliances 
only. 

The  necessity  of  transient  temperature  tests,  in  order  to 
isolate  the  exact  vertical  shift,  and  thus  obtain  maximum 
predictive  power,  should  be  explored. 

Reversibility  is  always  dependent  on  the  choice  of  a  well- 
defined  easily  reproducible  reference  state.  In  composites, 
however,  initial  stresses  arise  as  a  result  of  physio- 
chemical  processes,  technological  operations  and  mechanical 
and  thermal  incompatibility.  The  natural  state,  i.e.  a  body 
free  of  stresses  and  strains,  is  nonexistent.  The  importance 
of  an  "eigenstress"-field  on  the  viscoelastic  behavior  should 
be  checked. 
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APPENDIX  A 


Cure  Procedure  for  T300/934  Composite  Panels 


1.  Vacuum  bag  -  insert  lay-up  into  autoclave  at  room  temperature. 

2.  Apply  full  vacuum.  (A  minimum  of  25  inches  of  Hg.)  Hold  for  30 
minutes  at  room  temperature. 

3.  Maintain  full  vacuum  throughout  entire  cycle. 

4.  Raise  temperature  to  250  F  (+5 .  -  10  F)  at  2  ±  5  F  per  minute. 

*5.  Hold  at  250  F  (+5.  -  10  F)  for  15  ±  5  minutes,  apply  100  psi 

(+5.  -  0  psi ) . 

*6.  Hold  at  250  F  (+5.  -  10  F)  and  100  psi  (+5.  -  0  psi)  for  45  ±  5 
minutes. 

*7.  Raise  temperature  to  350  F  (+10.  -  0  F)  at  2  ±  5  F  per  minute. 

*8.  Hold  at  350  F  (+10.  -  0  F)  for  two  hours  ±  15  minutes. 

*9.  Cool  under  pressure  and  vacuum  to  below  175  F. 


*Pressure  was  85  psi . 
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Curing  Procedure  for  Resin  Panel  #1 


1.  Resin  heated  until  fluid  and  poured  into  heated  mold. 

2.  Mold  vacuum  bagged  and  inserted  into  autoclave. 

3.  Full  vacuum  applied  and  held  30  min. 

4.  Raise  temperature  to  250  F  (+5.  -10  F)  at  2  ±  5  F  per  minute. 

5.  Hold  at  250  F  (+5  -10  F)  for  15  ±  5  minutes.  Apply  100  psi 

(+5  -  0) 

6.  Hold  at  250  F  (+5  -10  F)  and  100  psi  (+5  -  0  psi)  for  45  ±  5  min 

7.  Raise  temperatuer  to  350  F  (+10  -0  F)  at  2  ±  5  F  per  minute. 

8.  Hold  at  350  F  (+10  -0  F)  for  two  hours  ±  15  minutes. 

9.  Cool  under  pressure  (85  psi)  and  vacuum  to  below  175  F. 


Fabrication  Procedure  for  Resin  Panel  #2 


1.  Mold  heated  to  200  F. 

2.  Resin  heated  to  160  F  and  degassed  in  bell  jar  for  5  minutes. 

3.  Degassed  resin  poured  into  one  corner  of  tilted  and  heated  mold. 
Air  gun  used  to  keep  mold  warm. 

4.  Resin  then  cured  in  oven  with  same  schedule  as  prepreg  tape,  i .e 
steps  4-9  above. 
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The  Curing  Procedure  for  Resin  Panel  #3 


934  Resin  Cure  Procedure  -  Solution  Form 


1.  Break  off  chunks  of  solid  (cold)  resin. 

2.  Place  in  beaker  of  acetone  or  methyl  ethyl  keytone. 

3.  Slowly  warm  and  stir  solution  up  to  80°F  (26°C)  and  until  all 
resin  is  dissolved  to  taffy  consistency. 

4.  Pour  into  teflon  coated  bread  pan  and  place  in  vacuum  oven. 

Heat  to  150°F  (66°C)  under  full  vacuum.  Observe  until  bubbles 
cease  (about  2  hours) . 

5.  Remove  from  vacuum  oven  and  put  in  212°F  (100°C)  oven  over  night 
(roughly  3:00  p.m.  -  7:00  a.m.). 

6.  Remove  from  oven  in  morning. 

7.  Put  in  oven  at  325°F  -  350°F  (163°C  -  177°C)  for  two  (2)  hours. 

8.  Turn  off  oven  and  allow  to  cool  gradually. 
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appendix  b 


APPENDIX  C 


INFLUENCE  OF  EXCENTRIC  LOAD  APPLICATION  [105] 

Suppose  that  the  load  F  is  applied  with  a  small  excentricity  e. 
This  causes  not  only  the  homogeneous  stress  state  which  we  intended 
to  produce,  but  also  the  bending  moment  F-e  as  is  shown  in  Fig.  C.l. 


Fig.  C.l  Superposition  of  Normal  Stress  and  Bending  Due  to  Excentric 
Load  Application 


The  stress,  at  the  gage  location  can  then  be  written  as, 

„  _  F  F-e 

aGL  w*t  ,2 
w  •  t 


=  a 


1  +  6* 


with  au  =  ^  ,  the  uniform  stress,  which  we  intended  to  apply. 
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